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SECTION 1 


INTRODUCTION AND SUMMARY 


I.I. Introduction 


This interim report is submitted in partial fulfillment of documen¬ 
tation requirements for the study: Active Control of Space Structures-Six 
(ACOSS Six). For the successful control of large space structures (LSS), it 
is recognized by the LSS community at large that the fundamental problem is 
the design of a finite-dimensional compensator to control an infinite 
dimensional system (1). To date, numerous theoretical contributions towards 
"solving" this problem have been made; however, the resolution of this 
problem requires further theoretical research which must be validated through 
appropriate design, analysis and experimentation. As a result of our endeavors, 
together with exposure to a broad spectrum of knowledge provided by the LSS 
community, it is our judgment that the technical issues in LSS control 
technology include the following: 

(1) LSS Modeling accuracy should be known to within some specified 
bounds. Modeling errors will limit achievable control system 
performance. These errors may be introduced through initially 
assumed structural properties or the truncation process implicit 
in the finite-element method. In space, LSS parameters may 
vary as a function of thermal gradients, configuration changes, 
or depletion of consumables. The more stringent the mission 
performance requirements, the greater the LSS model fidelity 
required. 

(2) Upper Atmosphere Models must be improved and verified by appro¬ 
priate experiments. Accurate knowledge of the external forces 
(e.g., earth magnetic and gravitational fields, solar wind and 
radiation pressure, drag) acting upon a LSS is necessary to 
satisfy precision control requirements. 

(3) System Identification is necessary for the purpose of LSS struc¬ 
tural model verification. Parameter identification techniques 
must be developed such that modal frequencies, damping ratios, 
and mode shapes can be accurately determined. Consideration 
must be given to the type of sensors, onboard r*.ocessing 
requirements, data reduction, post-processing requirements, and 
control law configuration. 

(4) Sensor and Actuator specifications must be determined in order 
to assess the applicability of existing hardware, as well as to 
provide new directions in research and development. 

(5) Control Law Design Methodology must address the following: 

(a) The model reduction process which reduces the high¬ 
dimensional finite-element model to a lower order design 
model. 
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(b) The design method for determining reduced-order 
compensators. 

(c) The criteria for determining overall closed-loop 
system stability. 

(d) Direct digital design methodologies and implementation 
techniques. 

(e) Sensor and actuator placement techniques which yield 
maximum observability and controllability. 

(f) Sensor and actuator dynamics. 

(g) The dynamic interaction of the attitude, vibration, and 
shape control laws. 

The theoretical research documented in the present report addresses control 
law design methodologies with particular emphasis on model reduction, design 
methods for determining reduced-order compensators, and criteria for 
determining overall closed-loop system stability. 

1.2 Research Scope 


1.2.1 Scope of Theoretical Research 


The enabling technology necessary for successful active control of 
large precision space structures requires mastery of at least the following 
principal topics: 

(1) Reduction of high-order, high-fidelity structural models to lower- 
order models suitable for the design of feedback controllers. 

(2) Design of fixed-order compensators which ensures stability for 
the total closed-loop system. 

(3) Analytical techniques for determining overall closed-loop 
system stability (including robustness). 

The theoretical results documented in this report address these topics. 
Given a finite-element modal representation of a LSS, a technique which 
aggregates the vibration modes with respect to actuator and sensor influences 
has been developed (Section 4). By proper selection of the aggregation 
transformation, this reduced-model may be made to approximate the open-loop 
performance of the LSS more accurately than methods utilizing direct modal 
truncation. New results relative to the extended Kosut design approach to 
multivariable output feedback controller design are presented in Section 2. 

The stochastic output feedback controller (Section 7) is formulated as the 
zero-order limit of the minimum variance fixed-order output feedback 
compensator. Section 5 discusses pareto-optimal control (multiple-objective 
optimization) with respect to the design of LSS vibration controllers. 
Stability margins for observer-based, fixed-order compensators is investigated 
in Section 3. The general conditions for LSS closed-loop stabilization by 
output feedback control are derived (see Section 6). 
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1.2.2 Scope of Applications Research 

In order to assess performance, sensitivity, and hardware requirements 
of the various active structural control methods, a simple but realistic 
evaluation model is necessary. Section 9 describes such a model (Draper 
Model //2) which is based on realistic sizes and masses. This model contains 
a simple optical system with unclassified performance measures and tolerances 
A line-of-sight (LOS) performance measure (Section 8) is available for 
evaluation purposes. Section 10 presents preliminary results of a sinusoidal 
disturbance (see Section 9) impact on the LOS performance. 

1.3 Summary of Principal Results 

1.3.1 Results from Theoretical Research 


Previously reported results relating to the algebraic consistency of, 
and the structure of solutions to, the Kosut gain equation are substantially 
improved (Section 2). Conditions for design-model-independence of free 
parameters occurring in such solutions are given. The nature and extent 
of performance degradation to be expected due to decentralizing the feedback 
information structure is obtained. Stable designs obtained by the extended 
Kosut method are shown to be robust. Major theoretical problems that need to 
be addressed are: (1) Characterization of systems for which the extended 
Kosut approach leads to a stable closed-loop design (at the design model 
level); (2) Synthesize a systematic procedure for assigning values to the 
free parameters (ut « >. available) so as to improve stability and/or perfor¬ 
mance of the full-order system incorporating the reduced-order controller. 

The stability margins for fixed-order compensators (Section 3) was 
investigated. The observer-based fixed-order compensator was assumed since 
it is representative of the most general case. The objective was to find 
a procedure which would determine the gain matrices such that the full-order, 
closed-loop system remained stable. If the resulting equations could be 
solved analytically, then explicit stability conditions could be given. Due 
to the complexity of these equations, an analytical closed-form expression is 
not apparent at this time. 

A standard approach to LSS model reduction is to truncate the open- 
loop modal representation directly by neglecting certain modes. The aggregate 
model development (Section 4) provides a technique for aggregating the 
vibration modes (state variables) with respect to actuator and sensor 
influences. The aggregation is accomplished by a similarity transformation 
which can be implemented by Gaussian elimination or Gram-Schmidt orthogon- 
alization procedures. By proper choice of the aggregation transformation, 
the resulting model may be reduced to a lower order while retaining an 
approximate open-loop performance equivalent to that of the full-order model. 
Since the aggregation transformation is nonunique, several approaches were 
examined which utilized this nonuniqueness for imposing additional properties 
on the aggregate design model. The feasibility of these approaches was also 
analyzed. 
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Pareto optimality (Section 5) is a practical notion of "multiple- 
objective" optimality that is useful ^uLrol designers faced with multiple 
disparate performance obieoM-. A design is Pareto optimal if none of the 

multiple objectives can be further improved without degrading any other. 

The education of Pareto-optimal control to the design of LSS vibration 
controllers provides the control designer with a systematic procedure for 
making efficient design tradeoffs. The underlying multiple-objective design 
problems relevant to LSS active control are shown via Canavin's modal dash- 
pot concept and Aubrun's low authority control theory. Favorable results 
were obtained by Pareto optimal designs on the 12-mode Draper tetrahedral 
model. Efficient tradeoffs between high damping and low output feedback 
gains were demonstrated. 

General conditions for closed-loop LSS stabilization (Section 6) by 
output feedback control were developed. If these conditions are satisfied by 
proper design, asymptotic stability of the closed-loop LSS system is ensured. 
The resulting theorems are applicable to the general case of LSS that contain 
not only stable elastic modes but also rigid-body modes and unstable modes. 

A specific output feedback control design technique is also presented for 
integrated control of elastic and rigid-body modes which increases both the 
damping and frequency of selected modes. This controller is closed-loop 
asymptotically stable as well as robust with respect to modal truncation and 
parameter variations. A technique is presented for synthesis of actuator in¬ 
fluences such that control spillover to selected modes is prevented. It is 
subsequently shown that this concept is readily extendable to placement of 
actuators which minimizes control spillover. 

The minimum variance fixed form output feedback compensator has been 
proposed as a possible solution to the problem of designing reduced-order 
controllers for plants of large dimensions. The zero-order limit of this 
compensator is the stochastic output feedback controller (Section 7). A 
detailed discussion of the stochastic output feedback controller and its 
relation to several other compensator formulations such as the minimum order 
dual observer based compensator and the Levine-Athans-Johnson reduced-order 
compensator is presented. Lower and upper bounds on the optimal stochastic 
output feedback controller cost are derived. These bounds should be useful 
in preliminary evaluations of subsequent controller performance. These 
bounds may also provide useful starting estimates in iteratively computing 
the optimal gains. Presently, this computation represents the major 
theoretical problem with the stochastic output feedback controller. 

1.3.2 Results from Applications Research 

The Draper model #2 (Sections 8, 9) was evaluated with respect to the 
sinusoidal disturbances given in Section 9. The preliminary results 
(Section 10) show that the addition of damping to individual modes is not 
sufficient to suppress a sinusoidal disturbance impact on the steady-state 
LOS amplitude error. 
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SECTION 2 


SUBOPTIMAL OUTPUT FEEDBACK - EXTENDED KOSUT METHOD 


2.1 Introduction 

2.1.1 Synopsis of Previous Research 

The Kosut approach to multivariable output feedback controller design 
was originally developed [1] as an alternative to optimal output feedback 
approaches [2] in which an acceptable solution to the strongly coupled 
nonlinear equations for the feedback gains was seldom obtainable [3]. By 
accepting suboptimal values for the control system performance measure, 
equations for the feedback gains were obtained that could be solved without 
iteration. Part of the price paid for the resultant simplicity of the 
design procedure was the lack of a guarantee for stability of the closed- 
loop system. As a result, the method found only limited application [4], 

Recently, a fresh look at the Kosut approach was taken to determine 
its feasibility for use as a design tool for controllers of large flexible 
spacecraft [5]. It was immediately recognized that certain assumptions on 
which the method is based are incompatible with corresponding properties 
of reduced-order modal models that are needed for controller design with 
large-scale systems; in particular, the state-to-output observation matrix 
in the reduced-order model must have maximum rank. It was observed that 
this restriction on the design technique is not in fact necessary. The 
linear equation for the feedback gain matrix can be solved in all cases; 
in particular, when this equation is rank-deficient, the general solution 
for the feedback gains contains free parameters proportional in number to 
the rank-deficiency. It was demonstrated how such free parameters can be 
selected so as to improve the performance of the full-order closed-loop 
system relative to what would be possible if the rank restriction on the 
observation matrix were retained [6]. As part of a parallel investigation, 
the Kosut approach incorporating the above extensions was compared to other 
representative multivariable controller design approaches in the context 
of a specific example having certain of the principal features of a large 
flexible spacecraft model [5,7]. In spite of the lack of a guarantee for 
stability of the closed-loop design model, the performance of the example 
system incorporating a Kosut controller compared very favorably with that 
using other controller designs. Finally, an algorithm for the numerical 
solution of rank-deficient linear equations was developed [8] to enable 
application of the extended Kosut approach to high-order design models. 

2.1.2 Summary of New Results 

New results relative to the extended Kosut design approach appearing 
in the present report are briefly summarized. A general criterion for 
stability of a closed loop system - at the design model level - for output 
feedback systems is given in terms of the basic parameters of the full-order 


6 







structural model, including sensor and actuator parameters, and the 
feedback gain matrix. A simple and useful, but undesirably restrictive, 
sufficient condition for stability is also given. Robustness of stability 
is established for designs obtained by the extended Kosut method in the 
situation that only rate sensing is used. In addition, certain previously 
reported results are improved. In particular, conditions are given under 
which free parameters occuring in the feedback gain matrix do not influence 
the stability of the design model. Under such conditions, a decision on 
their selection can be based solely on considerations of stability and 
performance of the full-order closed-loop system. Also, a simplified proof 
is given of the general algebraic consistency of the linear equation for 
the feedback gains [6]. 

It has been noted previously [5] that the essential nature of the 
Kosut design approach, including extensions, is not fundamentally altered 
when a decentralized information structure is imposed in the feedback path. 

A closer look at the differences in detail, however, suggests a much more 
satisfactory way of representing and understanding the general solution to 
the rank-deficient gain equation. A modified representation for the general 
solution is presented that is more readily related to the classical algebraic 
theory of linear spaces. The algorithm for numerical solution of such 
equations is correspondingly modified. A much clearer understanding of the 
impact of decentralization is revealed by a study of its implications upon 
performance for a specific system. Performance of the two-mass oscillator 
system incorporating several decentralized Kosut designs is examined. 

Specific limitations on system performance are observed when compared with 
corresponding centralized designs. 

The section is concluded with a brief statement of important questions 
that remain to be resolved regarding the extended Kosut approach. 

2.2 Notation 


In order to be self-contained, the notation to be used is briefly 
summarized. It is assumed that the structural dynamics of systems to be 
discussed can be adequately represented by a finite-element model: 

Mq + Kq * f (2-1) 


T 

where q = (q^,...,q n ) is the vector of generalized physical coordinates, 

M:n x n is the (positive definite) mass matrix, K:nxn is the (positive 

T 

semidefinite) stiffness matrix, and f = (f,,...,f ) is the vector of forces 

1 n 

applied to the structure. Equation (2-1) may be transformed to an equivalent 
representation in modal coordinates: 

ri + 0 2 n - <t T f = u (2-2) 

A 
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2 

where ft :n*n is the diagonal matrix containing squares of the natural 
2 2 

frequencies o)^,...,c4>:nxn is the nonsingular matrix whose columns are 
the structural mode shapes ... , 4> n , normalized such that 

$ T M4> =» I , 4> T K4> = ft 2 (2-3) 

n 

T 

B.rnxm is the actuator influence matrix for the control u = (u,,...,u ) , 

A T 1 m 

and n = (n 1 ,...,n ) is the vector of modal coordinates defined by 

q = 4m (2-4) 

Superscript "T" denotes matrix transpose, and the notation 1^ denotes the 
n*n identity matrix. The linear measurement equation has the form: 

y = c p <l + c v 9 (2-5) 

T 

where y = (yj,...,y^) is the vector of system outputs, and C R :£.xn, Cy:£xn 
are the displacement and rate measurement coefficient matrices, respectively. 

The structural model (2-2) may be placed in a convenient state space 
form as follows. Without loss of generality, it can be assumed that the 
components of n have been reordered so that the first N (N $ n) of them are 

T 

the selected critical modes n c = (r| »- - - » r » N ) » and that the remainder are 

residual modes n D = (i\,.... ) T . Corresponding submatrices ft 2 :NXN, 

ft R :(n-N)x(n-N) of u and 4> c :nXN, 4> R :nx(n-N) of 4> are defined similarly. 

With x = (n c , n c , t| R » n R ) as a state vector, the structural model (2-2), 
(2-4), (2-5) becomes: 

x = Ax + Bu (2-6) 

y = Cx (2-7) 
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Using residual mode truncation as the.model reduction process, 
order model with the state x c . (n c » Oq) is obtained: 


A_x_ + B„u 


a reduced' 
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New Results for Centralized Controllers 
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2.3.1 Stability and Robustness for Output Feedback Designs 

2.3.1.1 General Criteria for Stability 

As noted previously, one of the principal objections to the use of 
the extended Kosut method for controller design is the absence of any 
useful conditions which guarantee closed-loop stability even at the design 
model level. The approach taken to developing such conditions has been to 
try to discover an appropriate subclass of the systems (2-8), (2-9) being 
considered for which the extended Kosut approach leads to a stable closed- 
loop system at the design model level. 

The first result characterizes the eigenvalues of the closed-loop 
design model for general output feedback systems in terms of fundamental 
model parameters. 

Theorem 2-1 . Assume that the reduced-order model (2-8), (2-9) is connected 
by an output feedback law: u « Gy^,. ^ en eigenvalues of the closed- 

loop system are the roots of the (2N)th degree polynomial equation: 

det | ♦ C T [-A 2 M + *B A GC y - (K - V^p)] * c ( » 0 D (2-10) 


The proof is given in Section 2.5. A remark about Equation (2-10) is 
appropriate. Since the mode shapes $ ,...,$ n are linearly independent, 
the matrix = [ 0 ( 4> ] :n*N has maximum rank. It does not follow that 

equation (2-10) can be replaced with the (2n)th degree polynomial equation: 

-A 2 M + AB A GC y - (K - B A GCp) ] - 0 (2-11) 

The latter equation is easily recognized as the closed loop characteristic 
equation for the full -order structural model (2-1), (2-5) connected with the 
output feedback law u *= Gy. (Set A - ico in equation (2-11) to obtain an 
equivalent, more familiar, form.) Equation (2-11) has 2(n - N) more 
solutions than equation (2-10). 

A simple, but useful, sufficient condition for closed-loop design 
model stability is based on Theorem 2-1. For clarity, the result is 
presented for the case where rate sensors only are used (Cp * 0). 

Corresponding results for the general case are easily obtained. 
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Theorem 2-2 . Assume a closed-loop structure as In Theorem 2-1 with C * 

T * 

Suppose that there exist 6^ > 0, i - such that * diag (-6^). 

Then the closed-loop system is stable. Q 

The proof is given in Section 2.5. This result is useful for establishing 
robustness of stability. 


Specializing from general output feedback systems to those designed 
using the extended Kosut method, one may ask under what conditions the 
feedback gain matrix generated by that method satisfies the hypotheses of 
Theorem 2-2. This question has not been answered in general. It should 
be evident, however, that Theorem 2-2 is by no means necessary for design 
model stability; it is in fact a quite restrictive condition. This may be 
seen by examining the condition in the context of the two-mass oscillator 

T 2 T 2 

example worked out previously [5;Sec. 6]. In that case, <J> C ■ ($ ) G<J> , 

a scalar (hence diagonal), where 


G - G(a;e,6) 


*2 

-a + e /<J> 


-5 


♦ l 

-a + 6 V<f> 


The stability condition of Theorem 2-2 reduces here to the condition that 
G be negative definite, which requires [10] that: 


and 


y 2 2 

O - € /$. > 0 


[- 


y 2 2 y i 2 

e - 6 V*; 




However, it was seen [5;Sec. 6] that this closed-loop design model A. + B C GC C 

is independent of the free parameters e and 6, and is in fact stable as 
long as a > 0. 
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2.3.1.2 Robustness of Stability for the Extended Kosut Method 

The robustness question is examined in two parts. First, the 
sensitivity of the design model closed-loop eigenvalues is examined relative 
to changes in the feedback gain matrix. This is applicable to any output 
feedback system. Second, the sensitivity of the gains generated by the 
extended Kosut method is examined relative to changes in the fundamental 
structural model parameters. For clarity, it is again assumed that rate 
sensing only is used (Cp = 0). 

Theorem 2-3 . Assume a closed-loop structure as in Theorem 2-1 with C p = 0. 

The mapping which assigns to each output feedback gain matrix the set of 
closed-loop eigenvalues is continuous. D 

The proof is given in Section 2.5. The significance of Theorem 2-3 is 
as follows. Suppose that for a particular gain Gq the system is stable. 

Then the corresponding eigenvalues E (X*j*,... jX^) lie in the Cartesian 
product of the sets ReCA^) < 0, an open set in real Euclidean 2N-space. 

Theorem 2-3 says that there is an open set containing Gq in the space of 

m*£ real matrices for which the corresponding eigenvalues also satisfy 
Re(A^) < 0, i = 1,...,2N. It should be noted that Theorem 2-3 is 

essentially without significance if at least one stable set of eigenvalues 
cannot be found (cf. Theorem 2-2). 

For the following result, it is assumed that structural model reduction 
is by truncation (equations (2-8), (Z-9)), and that a specific reference 
system and a specific suboptimal performance measure for the extended Kosut 
method (minimum error excitation) have been chosen. The mapping described is 
then well defined. 

Theorem 2-4 . Assume C p = 0. The mapping which assigns to each mass¬ 
stiffness matrix pair in Equation (2-1) the extended-Kosut feedback gain 
matrix (Moore-Penrose inverse particular solution [8]) is continuous. 0 

The proof is given in Section 2.5. Together with Theorem 2-3, this result 
shows that if^ a stable design can be achieved by the extended Kosut method, 
then stability is retained under appropriately small changes in the 
fundamental structural model parameters. 

2.3.2 Refinement of Previously Reported Results 

Extensions to the Kosut method reported previously [5] are based on 
the existence and properties of solutions G:mx£ to the Kosut gain equation: 

g(c c PC c T ) - F*PC c T (2-12) 
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where F*:m*2N is the state feedback gain for the stable reference system used 
in the suboptimal design, and P:2Nx2N is the positive definite Lagrange 
multiplier matrix obtained from the Liapunov equation: 

(A c + B C F*)P + P(A C + B c F*) T + I 2n = 0 (2-13) 

It was shown [5;para. 6.2.3.3] that Equation (2-12) is solvable regardless 

T 

of the rank of (equivalently, rank of C^PC^ ), and that in the rank- 
deficient case (say, rank = r < £ £ 2N), the general solution has the form: 

G(D = G q + TS (2-14) 

where Gq is a particular solution of Equation (2-12), T is a matrix whose 

m*(i.-r) non-zero elements are free parameters, and S is an appropriate non¬ 
singular matrix. 

In a detailed design example, it was observed that the eigenvalues of 
the closed-loop design model were independent of the (two) free parameters 
in the feedback gain matrix. Since the performance (as contrasted with the 
stability) of the design model is not of practical interest, it was possible 
to base the selection of free parameter values solely on full-order system 
performance considerations. The question arises as to what conditions 
characterize this "independence" in general. To begin, a very simple 
sufficient condition can be stated. 

Fact 2-5 . Assume a closed-loop structure as in Theorem 2-1; refer to 
Equation (2-14). If B^TSCq = 0, then the eigenvalues of the closed-loop 

design model are independent of the free parameter matrix T. 0 

The proof is self-evident, since both G(T) and Gq generate the same closed- 

loop system matrix. It is worth observing that this condition (in fact, 
even the stronger condition: TSCq « 0) has been satisfied in all numerical 

examples studied to date. The following theorem fully characterizes the 
independence condition. 

Theorem 2-6 . Assume a closed-loop structure as in Theorem 2-1; refer to 
Equation (2-14). The eigenvalues of the closed-loop design model are 
independent of the free parameter matrix T if and only if Aq + B C G(DC C and 

A c + BqGqCq are similar; i.e., if and only if there exists a nonsingular 
matrix Q:2N*2N such that: 


[A c + B c G(r)c c ]Q - Q [A c + B c G o c c ] * 0 • D 
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The proof is given in Section 2.5. The full implications of this freedom 
of choice in selecting a feedback gain matrix with the extended Kosut 
method have not yet been uncovered. Performance improvements for low-order 
vibrating systems have already been demonstrated [5,6]. It is worth 
observing, however, that the question of effective use of available freedom 
of choice in selections of feedback gains is an important, and largely 
unexplored, question for general multivariable feedback systems [13;p.51]. 

Proof of the general algebraic consistency of Equation (2-12) has 
been considerably simplified (cp. [6;Theorem 2], [5;Theorem 6-6]). Since 
this is a principal result, the streamlined proof is of interest. Details 
are given in Section 2.5. 

2.4 New Results for Decentralized Controllers 
2.4.1 New Perceptions Suggested by Decentralization 


The discussion of decentralization in this section of the report is 
restricted to the nature of the information structure in the sensor-to- 
actuator feedback connection. Totally decentralized systems which include 
decentralized plant structures are not treated. Specifically, it is assumed 
that the information for driving a given actuator is derived from a (possibly 
proper) subset of all available sensors, and that the collection of sensors in 
this subset is (in general) different from one actuator to another. In 
symbols, the sensor combinations for transmission to each actuator are 
generated from the state vector as follows: 

y 1 - C^x , C i :«, 1 x4N 5 i = 1.m (2-15) 


Equation (2-15) replaces Equation (2-7) in a decentralized full-order model. 
The feedback connection is made by: 


T i 

= g*y 


VV 1 


i = 1,...,m 


(2-16) 


Remark : In the special case that the components of each y* are selections from 

_ t 

the components of a fixed vector y =(y^,...,y^) , Equation (2-16) may be 

_ T 

written in the form u = Gy, where u = (u,,...,u ) , and G = [g..] has the 

l m ij 

property that g = 0 if the component y of the fixed vector y is not also 
ij i J 

a component of the output vector y associated with u^. In general, however, 
this situation will not necessarily exist. 
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The Kosut gain equations corresponding to Equation (2-12) for a 
controller with decentralized feedback structure are [5]: 


T 

8 i 


( c c lpc J 


) 


•ft i 

F* p C c 


i = 1,... ,m 


(2-17) 


* * i 

where F^:lx2N are the rows of F , and C^:Jl^x2N are the matrices corresponding 

to Equation (2-15) for a reduced-order model. The structure of solutions to 
algebraically consistent equations of the form (2-17) may be expressed in 
the following form. 

Theorem 2-7 . Let A:A*V, b:lxv be matrices for which the equation xA = b 
has a solution x^. Denote r = rank (A). The general solution of this 

equation has the form: 


x 


T 

+ y S 


(2-18) 


where the A-r row vectors of S:(A-r)xA form a basis for the null space of 
A (viewed as a mapping of E^ into E V ), and the (A-r)-vector y is arbitrary. Q 

The proof is given in Section 2.5. This result immediately suggests an 
appropriate reformulation of previously reported solution structures for 
homogeneous matrix equations [cf. 5, Theorem 6-3; 6, Theorem 3; 8, Theorem 3-9] 
which sharpens the statement and identifies the true nature of the "S" 
matrix. 

Theorem 2-8 . Let A:Axv be a matrix. Denote r = rank (A). The general 
solution of the equation XA = 0^ has the form X = TS, where the A-r row 

vectors of S:(A-r)xA form a basis for the null space of A (viewed as a 
mapping from E^ into E V ), and the matrix r:yx(A-r) is arbitrary. D 

This result shows that the last A-r rows of the S matrix in Theorem 6-3 of 

T 

[5] and the last Z-r rows of the U matrix in Theorem 3-9 of [8] correspond 
simply to different choices of a basis for the null space of A. 

Step 4 of the algorithm for numerical solution of the Kosut gain 
equation [8, Sec. 3.5] should now be modified so as to be consistent with 
Theorem 2-8. 

Step 4 (revised) . Find the general solution. 

(a) Select an m*(£-r) matrix f of free parameters. 
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(b) Choose row Jl-vectors s^,. 
space of A and form the matrix 


. ,s„ that form a basis for the null 
£,-r 


S & 


Z- r 


(One possible choice is the set of vectors obtained by transposing 
columns r+l,...,£ of the matrix U from Step 1(a)). 

(c) Compute the homogeneous solution X = TS. 

c 


(d) Compute the general solution X = X + X 
Steps 3(b) and 4(c). Q P c 


using the results of 


Step 4(b) makes explicit the fact that there is considerable freedom in 
choosing the matrix S (although not complete arbitrariness). Specific 
approaches to exploiting this freedom of choice await further investigation. 


2.4.2 Performance Implications of Decentralization 

There are a number of reasons why a decentralized information structure 
might be incorporated in a physical system. For a very large structure with 
many sensors, the computational burden of incorporating information from 
all of the sensors in the determination of control signals to each actuator 
may far outweigh the advantages when a centralized processing unit is used. 

If numerous processing units are distributed throughout the structure, it 
may not be feasible to transmit information from all the sensors to each 
processing unit. On the other hand, there are cases in which the only 
relevant information is that from nearby sensors (e.g., in the case of 
shape control), in which case information from all the sensors at each 
processing unit would be largely unnecessary. Whatever the reasons, however, 
it is important to take these information structure constraints explicitly 
into account in the mathematical model used for controller design and eval¬ 
uation. In general, one would expect some loss in performance (by whatever 
measure) as the price of operating with restrictions on the internal 
information structure for determining control signals. 


A clear picture of such performance implications does not emerge from 
the design equations themselves. In fact, the gain equations (2-17) for the 
decentralized controller have exactly the same mathematical structure as the 
gain equation (2-12) for the centralized controller. Hence, all the results 
obtained for the latter (e.g., algebraic consistency, structure of general 
solution) apply to the former. Therefore, some design calculations were 
made on a specific representative system (the two-mass oscillator example [5]) 
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to get at least a snapshot of the nature of performance degradation, if any, 
associated with decentralization of the controller. The oscillator was 
assumed to have one force actuator and one rate sensor associated with each 
mass [5; Fig. 2-2]. Each of three possible types of decentralization in the 
feedback structure was studied: (a) fully decentralized: = u^(y^), 

u 2 = u 2^ y 2^* (k) partially decentralized with actuator 1 localized: 

U 1 = u l^ y l^’ u 2 = u 2^ y l’ y 2^ ; ^ P artiall y decentralized with actuator 2 
localized: u^ = u^y^, y 2 ), u 2 = u 2 <y 2 >. The reference system for the 

Kosut design was taken as determined in previous examples [5], Hence F*, 
and also P, [cf. Equation (2-13)] are as previously computed: 




0.066389774 

0.51609641 

I 

N3 * ' 

I 

i _: 


-0.028241581 

L 

-0.21954252 



1.1373892 

-0.5 


i 

= 

-0.5 

7.4375655 _ 

• 


2.4.2.1 Fully Decentralized Design 

Recall that the state vector for the reduced-order model 
[cf. Equation (2-8)] is x c = (f] 2> r^)* that residual mode (n^) dynamics are 

being ignored, and that Equation (2-4) relates the physical and modal 
coordinates. The specifications for this configuration may be stated: 

1*2* i 

y = q^, y = q 2> Accordingly, the matrices corresponding to 

Equation (2-15) for the reduced-order system are: 

C J ■ [° *l] • C C ' [° *2 2 ] • 

Since each C* has maximum rank, each of the gain equations (2-17) has a 
unique solution. The solutions are: 



where = 0.59708155 is the reference-system-dependent parameter occurring 
in the solution of Equation (2-12) for the centralized design [5, Eq.(6-43)]. 
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Referring to the remark following Equation (2-16), we may write the 
feedback connection as u = Gy^,, where 



The closed loop design model has the same characteristic polynomial as in 
the centralized case: 

X 2 + a°(<J ) 2 <p 2 )\ + ca 2 

It is stable, with u) = co^ (open-loop critical mode frequency), and £ = 0.1 

(design objective). Note that, in contrast to the corresponding centralized 
design, there are no free parameters in the solution for the feedback gains. 

The matrices C 1 of Equation (2-15) for the full-order system: 

c 1 = [o <J> 2 0 4)}], C 2 = [o <p 2 0 4>2^J 

are used to connect the reduced-order controller to the full-order system 
via Equation (2-16). The characteristic polynomial for the closed-loop 
full-order system is 



Dynamic characteristics of this system are: 

Mode 1 (critical): u)^ = 2.5851902 sec 

= 0.099939961 

Mode 1 (residual): = 0.54704429 sec 

he 
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0.34631083 





Some effects of decentralization on overall system performance can be seen 
by comparing the time response in Figure 2-1 and the frequency response in 
Figure 2-2 with the corresponding plots for the centralized design 
[5; Figs. 6-7 and 6-8, respectively]. Noticeable degradation does occur, and 
is traceable to the reduction in damping of the residual mode. Nonetheless, 
the degradation is not nearly as severe as that occurring when the sensor 
on the outer mass of the two-mass system is deactivated while retaining a 
centralized information structure on the remaining available sensor 
[5; Figs. 6-13 and 6-14, respectively]. 

2.4.2.2 Partially Decentralized Design: Actuator 1 Localized 

The specifications for this configuration may be stated: y^ = q., 

2 • • T i * 

y = (q^, ^2^ ' The re ^ atec ^ matrices C^, corresponding to Equation (2-15) 

for the reduced-order system are: 


r 21 

< 

<N 

r« , 2*1 

L° *ij * 

in 

1 u 
u 

L° * l 


Solutions to the gain equations (2-17) are as 


g x = -o° , g^(6) = |^-S 

where 5 is a free parameter. An equivalent feedback connection of form 
u = Gy^, exists and is: 



T 


0 


8 1 ° 


-a 0 

G(<$) £ 


= 

2 


g 2 (< 5) 


-6 -o° +6 1 !$\ 


follows: 

-a° + 6^ L /4>2 ] 


The characteristic polynomial for the closed-loop design model is the same 
as that for the preceding design; in particular, it is independent of the 
free parameter 6. 

The matrices C 1 of Equation (2-15) for the full-order system are: 


1 A 

r 

2 


1 "1 

2 A 

r 

, 2 


1 ~1 

C = 

t° 

♦l 

0 

M • 

C = 

L° 

<p 

0 

♦ J 
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To assist in determining an appropriate value for the free parameter 
the modal closed-loop equations are developed: 


l/ 2 T .2\* 2 

r o/. 2 T ,i\ r 

\4> 4> /n 2 + w 2 n 2 = 

1-0 (4) 4> / - 6 det 4> 1 i 

(J) 1 ) + 6 det 4> () 

l] rij + cajrij = [ ~ct°(<|> 2 


( 2 - 20 ) 


Several possibilities for selecting 6 are as follows: (a) choose 6 so that 
the coefficient of in Equation (2-19) vanishes; (b) choose 6 so that the 

coefficient of in Equation (2-20) attains some desired value; or 

(c) some combination of (a) and (b). Selection (a) eliminates the effect on 
the critical mode of residual mode dynamics, whereas selection (b) approx¬ 
imates a specification for the residual mode damping. Past design examples 
have indicated that failure to damp the residual mode strongly and adversely 
affects the full system performance, and that the spillover effect admitted 
by Equation (2-19) is a lesser concern. Since selection (a) will not affect 
residual mode damping, whereas selection (b) will tend to reduce the 
spillover problem indicated by Equation (2-19), selection (b) is made. Thus: 


P 0 A 


2c id"i ' “°(* 1 + 1 ) 

(+!/♦!) 


0.33580542 


det $ 


where is chosen to be 1 /-/l. The gain values then become: 


G(<5°) = 

. - 0 . 


59708155 

33580542 


-1.3864866 


With this connection, the dynamic characteristics of the system are: 


Mode 2 

(critical); 

“2C 

2.5914658 sec - 



C 2C 

0.10021320 

Mode 1 

(residual): 

W 1C 

0.54571957 sec' 



? 1C 

0.70633967 
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Performance of this design is so close to that of the corresponding 
centralized design that no degradation can be observed from the plots. 
Figures 2-3 and 2-4 may be compared with the corresponding plots for the 
centralized system [5; Figs. 6-7 and 6-8, respectively]. Consequently, 
Figures 2-3 and 2-4 may themselves be viewed as representing the 
centralized system time and frequency response, respectively, for the 
purposes of the comparisons in this section. 


Decentralized Desif 


Actuator 2 Localized 


The specifications for this configuration may be stated: 

1 . . t 2 i 

y = (q^t q 2 ) » y = q 2 - The related matrices corresponding to 

Equation (2-15) for the reduced-order system are: 


4 - [° ♦*]• 4 - [° ♦!] 


Solutions to the gain equation (2-17) are as follows: 


r a 2 

0 ^2 2 

-a u + e -ej , 


where e is a free parameter. An equivalent feedback connection of form 
u = Gy c is: 



The characteristic polynomial for the closed-loop design model is the same 
as for the fully decentralized design; it is also independent of the free 
parameter e. 

The matrices C 1 of Equation (2-15) for the full-order system are: 

c 1 = ^0 <J > 2 0 4) 1 ] , c 2 = [o 4> 2 0 4>2 ] 
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fQTEM: Cigga RESPONSE; Oa/Fg time- 17 • S 7:88.3 DOTE- D9/C9/80 180/283 



Figure 2-4. Frequency response of q£ to periodic disturbance f 2 (t) - sin wt 
(partially decentralized controller: actuator 1 localized). 




The modal closed-loop equations (to help select e) are: 


n 2 + o°(<f> 2 4> 2 )n 2 + ^2 * M* 2 + e det *] n L 

) i + [o 0 ^ 1 4* 1 ) - e det $ ( 1 /4> 2 )]n 1 + = [-o°(^ 2 (j) 1 )] n 2 


(2-21 


Following the same arguments as In the preceding section, the parameter e 
is chosen so that the coefficient of n, in Equation (2-21) has the desired 
value = 1//5: 


0 

e 


A 


2?^ - a 0 ^ 1 4* 1 ) 
det 4> (4>|/-<j) 2 ) 


0.92785918 


The gain values then become: 


G(e°) 




-0.99178405 

0 


-0.92785918 

-0.59708155 


With this connection, the dynamic characteristics of the system are: 


Mode 2 

(critical): 

u)_ ^ = 
2C 

2.5683141 sec" 



C 2C 

0.098354234 

Mode 1 

(residual): 

U 1C = 

0.55063889 sec' 



C 1C 

0.71291349 


The only performance degradation relative to the centralized design 
observable from the plots is in the frequency domain, and is quite minor. 
The frequency response shows slightly higher magnitudes near a) = u> 2 

(open-loop critical mode frequency), and slightly higher phase shifts for 
w 2 Figures 2-5 and 2-6 may be compared with the corresponding plots 

for the centralized system [5; Figs. 6-7 and 6-8, respectively]. 
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2.4.2.4 Remarks on Performance Comparisons 

A summary of the design studies on the two-mass oscillator example 
for all of the centralized and decentralized designs discussed in this 
section is shown in Table 2-1. Data for the centralized designs was 
previously reported [5,7]. 

Broad conclusions regarding the effects on performance of decen¬ 
tralizing the feedback information structure cannot be drawn from these 
few designs on such a low-order system. However, several observations are 
worth noting. First, noticeable performance degradation did occur, and the 
degree of degradation increased as the degree of decentralization increased. 
Second, the worst performance degradation (for the fully decentralized 
configuration) was not as bad as one might have expected. Finally, the 
design flexibility due to the free parameters occurring in the designs with 
partial, but less than total, decentralization was sufficient to eliminate 
(nearly) all performance degradation. 

2.5 Proofs of Theoretical Results 


The proof of Theorem 2-1 depends on the following two facts. 

USS-l- * "c 

Proof . It follows from a block matrix representation of Equation (2-3) 

that («{» i ) T M<J)^ = <5^, and ((p i ) T K<j>' ] * , *>3 “ 1. . where 6^ is the 

Kronecker indicator function: 

1 , 

0 , 

In particular, these relations hold for 

Fact 2 . (Rosenbrock [9]). If A, B, C, 
square and A nonsingular, then: 

T A 1 B 1 -i 

det-1-- det A det (-CA B + D) 

LC i Dj 



i = j 
i * j 

i, j = 1.N. I 

D are matrices, with A, D being 
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Proof of Theorem 2-1 . The closed-loop system matrix for the model (2-8), 
(2-9) under the connection u = Gy^ is: 


A c + B C GC C 


- a c + '!wc 


*5>a“»*c 


It is easily seen from the properties of determinants that the vanishing of 

2 T 

det (A^ + B C GC C ) is equivalent to the condition: det [-0^ + ^B^Cp^] = 0. 

This proves result (2-10) for the case X = 0. For X ^ 0, the matrix -Xl^ 
is nonsingular. Using Fact 2 and the properties of determinants leads 
to: 

det <A C + B C GC C - Alg,) ' d “ I ‘ >2l N + 4 C [B A G(XC V + V^C ' «cl 


Result (2-10) follows by an application of Fact 1. R 

Proof of Theorem 2-2 . Denote A = diag (-6^). Using Theorem 2-1 and Fact 1, 
the eigenvalues of the closed-loop system are **een to be the roots of: 


0 - det (-aVaa-O?) - <-» N n > 2 + V + “ 2 i) 


Since each 6^ > 0, all eigenvalues have negative real parts. 


Proof of Theorem 2-3 . Denote by G (- 


-> (X 1 ,...,y 2n ) the mapping 


described in the theorem. The function y is a composition ° ^ ° 

Yj t 

of the following three mappings: (1) G I — -» 

Y, 


(2) Z I—P z (*)> where P Z (X) = det ^X 2 ]^ + XZ + fi 2 \ ; and 
Y, 


(3) p(‘) h 


-> (X....X^^), the roots of p(X) = 0. It can be shown that 


Yp y 2> and y^ are continuous; hence y is continuous. Continuity of y^ is 

based on the well-known root-coefficient relations for polynomial 
equations [11]. I 
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Proof of Theorem 2-4 . Denote by (M,K) |-> G the mapping described in 

the theorem. The function 0 is a composition MM 0 of the following 

ft J ^ 1 

P 1 2 2 / 2 2 \ 
three mappings: (1) (M,K) | ——> (ft , $) , where ft = diag I oj ^.,u> ), 

uk being the roots of det (-w^M + K) =0, and $ = ^ j ••• j ^ n ] , <\> being 

the solutions of ^-u)^M + ^j(j) = 0 such that Equation (2-3) holds; 

2 ^2 

(2) (ft , $) |——> (A , B , C ), the latter matrices being defined by 

^ ^ ^ 3 

Equations (2-8), (2-9); and (3) (A c> B c> C c ) I——> G = (F*PC C T )(C C PC C T ) + , 

where F* is the state feedback gain for the Kosut reference system, and P 
is the solution of Equation (2-13). The superscript "+" denotes the 
Moore-Penrose inverse [8]. It can be shown that 0^, 0^, and 0^ are 
continuous; hence 0 is continuous. I 


The proof of Theorem 2-6 is self-evident from the following general 

result. 


Fact 3 . Two square matrices are similar if and only if they have the same 
set of eigenvalues. 

Proof . (Only if) This is a well-known result [12; p.144]. 

(If) Since both matrices have the same set of eigenvalues, they are each 
similar to the same Jordan block matrix. Since similarity is an equi¬ 
valence relation, the matrices in question are similar to each other. 


Simplified proof of Theorem 6-6 of [5; Sec. 6] , 


Denote E 


Since the rank of a matrix is equal to the number of 


T T 

its linearly independent rows: rank (C PC ) 5 rank (EC ). Conversely, 

suppose that this inequality is strict. Denote r = rank (EC r ). Using 

T ^ 

the fact that rank (C C PC^, ) = rank (C^) when P is positive definite, it 

follows that rank (C r ) < r. But then it follows easily [12] that 
T ^ T 

rank (EC ) - min | rank (E), rank (C p A ) } < r, which contradicts the 
definition of r. Hence, rank (C^PC^ ) - rank (EC^, ). I 


The proof of Theorem 2-7 is based on the following well-known result 
from linear algebra [12], 
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Fact 4 . Let <(■/ be a linear transformation defined on a vector space X. 
Then: 


dim [Range GV) ] + dim [Null CV) ] = dim (X) 


xA on the row-space e\ 


Proof of Theorem 2-7 . Define the map x 
Since rank ( t r/) = rank (A) = r, it follows from Fact 4 that Null (.V), 
null space of <V, has dimension A-r. If r = A, then Null (if/) = 0 and 
is unique. Otherwise, r < A, so there exist row vectors s 1t ...,s. 


A T T 1 | t T T 

that form a basis for Null (,«/). Denote S = s, i ••• i s, . 

i L 1 « i A-rJ 


r • ‘ ‘ 

T 

A- 


’A-r 


the 

X 0 


Let 


x e E satisfy xA = b. 


Then (x - x^) e Null (,V) , so there exist constants 


Y = (Y 


r A- 


) such that 


. A—r 

> Y.s. = 

^i=l 1 1 


T 

= y i s 


Conversely, any vector of form (2-18) satisfies the equation xA = b. I 
2.6 Conclusion 


A capsule summary of the results reported in this section may be 
stated as follows. First, certain previously reported results relating to 
the algebraic consistency of, and the structure of solutions to, the Kosut 
gain equation are substantially improved. Moreover, conditions for design- 
model- independence of free parameters occurring in such solutions are given. 
Second, a clearer idea of the nature and extent of performance degradation 
to be expected due to decentralizing the feedback information structure is 
obtained. Finally, stable designs obtained by the extended Kosut method 
are shown to be robust. 

However, some substantial open questions, both in theory and in 
application methods, must be resolved before the extended Kosut approach 
can become a competitive method for large space structure controller design. 

Major theoretical problems : 

1. Characterize systems for which the extended Kosut approach leads 
to a stable closed-loop design (at the design model level). 

2. Synthesize a systematic procedure for assigning values to free 
parameters (when available) so as to improve performance of the full-order 
system incorporating the reduced-order controller. 
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Major applications problems 

1. Develop the algorithm for solving the Kosut gain equation [8] 
into an effective design tool. 

2. Compare the extended Kosut method with several other represen¬ 
tative design methods on more sophisticated structural models. 
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SECTION 3 


STABILITY MARGINS FOR FIXED-ORDER COMPENSATORS 


3.0 Abstract 


In order that a theory can be developed for establishing an upper 
bound on closed-loop stability when evaluating a fixed-order compensator 
with higher-order evaluation models of large space structures, the closed- 
loop dynamic equations were derived first (see Subsection 3.1). As an 
initial effort, the common observer-based fixed-order controllers were 
considered; the design of such compensators was assumed to be based only 
on a reduced-order model formed by truncating all but a fixed number of 
vibration modes. The resultant closed-loop dynamic equations were then put 
into the standard second-order matrix form in two different ways, all having 
a certain form of effective damping and stiffness matrices. 

The closed-loop effective stiffness matrices were not symmetric; 
their symmetrization was then studied, so that standard stability conditions 
on second-order dynamics could be applied. Two steps were proposed: 
symmetrizing transformation of the coordinates and orthonormalizing 
transformation of the new eigenvectors. (See Subsection 3.2.) The closed- 
loop dynamic equations were then put in a modal form similar to the open- 
loop case. Determining the closed-loop stability and its margin would 
then be a straightforward matter after carrying out this two-step procedure, 
say, by a computer. Analytic expressions for the closed-loop stability 
conditions (with higher-order evaluation models), however, were not easy 
to derive. 

The stability of the closed-loop system having a fixed-order compen¬ 
sator could be looked at directly in the frequency domain. In other words, 
the shortest distance of the closed-loop poles to the imaginary axis of the 
complex plane could determine the upper bound of closed-loop stability. 

An example was worked out analytically in detail (see Subsection 3.3) to 
demonstrate that it would be very difficult to obtain an analytical 
expression of closed-loop stability margins from this viewpoint even for 
a structure as simple as having only two elastic modes. 
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3.1 


Closed-Loop Dynamic 


Compensators 


Equations with Observer-Based Fixed-Order 


Assume that an open-loop modal representation of the large space struc¬ 
tures, in the form 


i n + (G + D)n + ft 2 n = g u 
y = Pn + Vn (3-1) 

is available, and that the critical modes are the low-frequency ones while the 
residual modes are of high frequencies. G and D are gyroscopic and damping 
matrices respectively; ft is a diagonal matrix of natural frequencies. Let i, 

m, n, n , and n denote the numbers of outputs, inputs, modes, critical modes, 
c r 

and residual modes, respectively. Consider the following partitioning. 


“ ■ 
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n 
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, ft 2 = 

ft 0 
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cc 

cr 


cc 
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, D = D = 

d t 
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cr 

rr 
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[P P ), V 
c r 


[V. VJ 


(3-3) 


We presume that the parameters of the full model are available, but that the 
design is to be based essentially on the critical-mode part of the model. The 
feedback controls that will be considered are of the form 


u 


Kn + Ln 


(3-4) 
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where 


n +(G +D )n +fin 
c cc cc c c c 


8 u + H(y - P n - V n ) 
c c c c c 


In (3-5), n has the form of an observer for n . Note that it is driven by 
c c 

the output, y, which contains n r » and is thus not actually an observer of r> c 

in the formal sense. Evidently, (3-5) involves only the critical mode param¬ 
eters. The free design parameters are the gain matrices K, L, and H, and we 
want to know if they can be chosen by a procedure so that (3-1), (3-4), (3-5) 
yield a stable closed-loop system, particularly when G and D are very small. 

The reason for assuming this form is that it is representative of a 
range of practical schemes for finding gains. Note, for instance, that K, L, 
could be designed on the basis of an LQ control problem, and that H could be 
designed on the basis of a Kalman filtering problem or as a model matched 
filter for the critical modes. However, these methods have not been shown to 
lead in general to closed-loop stability, and more refined conditions on 
their application are needed. 

The standard coordinates for the closed-loop dynamics are n and 


e = \ \ (3-6) 

We write these in second-order form rather than state-space form. Let 



Effective damping 


38 






39 





3.2 Symmetrization of Closed-Loop Stiffness Matrices 

We note that the closed-loop system can be put into second-order form, 
but it is not clear whether this could be transformed to an equivalent de¬ 
coupled system 


+ [G + D] x + a 


2 ~ 


= 0 


(3-11) 


where 

-T ~ ~T ~2 

G = -G , D = D , and 0 is diagonal. 

If it could, then a rather simple stability test for second-order systems 

_2 

could be applied. For example [1], [2], when 0 is positive definite, system 
(3-11) is asymptotically stable if and only if the following associate system 

! z + Gz + fi^z = 0 

y = Dz 


is completely observable. 

The conversion of (3-10) to (3-11) proceeds in two steps: 

Step 1 : A symmetrizing transformation, S, which converts (3-10) 
so that its last term is symmetric. 

Step 2 : An orthonormal transformation, $, which puts the result of 
Step 1 into modal form. 

Step 1 

Suppose we perform a change of variable S as follows: 

Sv = x (S nonsingular) (3-12) 

For simplicity, rewrite system (3-10) in the form 

x + Fx + Kx = 0 (3-13) 

then 

v + S _1 FSv + S _1 KSv - 0 (3-14) 

If the last term is to be symmetric, there must be an S such that 
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-1- T-T -IT 

S KS = STS 


T-T T -1 

i.e., K = (SS )K (SS ) 


or equivalently, there must be a nonsingular symmetric transformation T such 
that 


KT = TK T 


In other words, such a transformation must be solved from the following mat¬ 
rix equation: 


“ HP T 11 T 12 


2 T 

-6K n t 12 t 22 


-2 T T 

T u t 12 a -KB 


T T T 2 

T 12 T 22 - PH 0 


ii 2 n 2 - B K + HP 
c c c 

To complete Step 1, suppose that such a symmetric nonsingular transformation 
T has been obtained. Then there is a nonsingular square root of T, denoted 
S, such that 


SS = T 


This can be found via Cholesky decomposition. 


K = S'^S 



Also define 





» 1 = »T 

G = -|(A - A 1 ) 

(3-21) 


= 1 - «T 

D = 2< a + A > 

(3-22) 

Then there exists an 

orthogonal transformation 4> such that 



T= -2 

4> K4> = ,1 

(3-23) 

Defining 

G = l> T Gi> 

(3-24) 


D = 4> T G0 

(3-25) 


brings the system into the form (3-11). 

Comment 1 : If steps 1 and 2 can be carried out, it is then possible 
to interpret the closed-loop system in terms of modes , i.e., as if it repre¬ 
sented a new structure. It was not previously apparent that "closed-loop 
modes" could be defined in the sense that open-loop ones were. 

Comment 2 : If any of the steps can be carried through analytically , 
then explicit stability conditions could be given. This has not yet been 
accomplished. 

Comment 3 : This procedure might represent a more simple way to 
determine closed-loop stability numerically using existing computer programs 
for high-order systems. 

Comment 4 : Since the design value of the foregoing results often de¬ 
pend crucially on the ability to obtain analytical expressions involving the 
gains, several examples were worked in order to obtain a better qualitative 
feel for how the gains enter the exact analytical expressions for stability. 
In symmetrization of a two-mode system with control of one critical mode, we 
carried out step 1 of the symmetrization procedures, Eq. (3-17), assuming 
n £ = n * 1. The expressions obtained were very complex. It should be ap¬ 
parent, however, that the problem of finding a symmetrizing transformation 
can be expected to yield interesting stability results. 

3.3 Closed-Loop Stability Conditions for a Two-Mode Example 

The stability problem of system (3-10) can be viewed directly in the 
frequency domain. In this case, we seek the roots of the polynomial 
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M(s) 


+ (C + D - t* I. + HV )s + ( 


K + HP ; 


det 


-H(Vs + F) 


-P(Ls + K) 


s + ((. + D)s + 


Note that the lower right block is just the polynomial of the open-loop 
system, which is stable by assumption. Thus a stability analysis might be 
based on the perturbation-theory expression 


M(s) = det { s 2 + (G + D - p L + HV )s 
| recce c 

+ (,. 2 - 8 K + HP ) 
c c c 

- H(Vs + P) [s 2 + (G + D)s + . 2 ] t( Ls + K) 
• det [s^ + (G + D)s + il ) 


An effort has been made to work out M(s) for a two-mode example, and to apply 
Routh's criterion so that exact stability conditions might be found for the 
simplest case of interest. The results follow 


M(s) 


fs 2 


+ (d 


6 

c 


S +hv )s+ 
c 


2 


- 6 k + hp ) 
c c 


det 


-6 (<s + k) 
c 

-6 r (S-s + k) 


-h(v s + p ) -h(v s + p ) 


s‘ + d b + (g + d )s 


2 2 

(-g + d ) s s + d s + ,. 


where lower-case letters have been used to denote scalar quantities. Note 
that g = g = 0 by the requirement that gyroscopic coupling be antisym¬ 
metric. Also d = d and g = -g have been used. Subscript c denotes 
rc cr °rc cr r 

"critical" and subscript r denotes "residual" as above. Lengthy algebra 
yields 


M(s) = £ M i si 

i-0 
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where 


-p i + v k + 2d +d 

c c cc rr 


h. (b v - B v ) + k(p + v Cd + d )) - £.0 (d + d ) 
r r c c c c cc rr c cc rr 

2 2 2 2 2 
-kt + g -d + 2d d +d + 2 cj + w 
c cr cr cc rr cc c r 


= hk(bv - S v ) + h!(B p +d B v -Bp - d Bv) 
rr cc rr ccrr c c rr c c 

2 2 
+h(d p + x v +d p +d d v + w v ) 
cc c c c rr c rr cc c r c 

+k(-d B - d B ) + t(-^ 2 B - d d B - u) 2 8 ) 

cc c rr c c c rr cc c r c 

2 2 2 2 
+2d _ + 2d . + d d + 2d ui 

cc c rr c rr cc cc r 


2 2 2 

hk B 6 p +hk6pp + hk(a B p + d Bv -Bp - d Bv 

err r r c c r r r cc r r c r c rr c c 

2 2 

+ht(d Bp + xu8v - d Bp - w B v ) 

cc r r err rr c c r c c 

+h[p (g 2 - d 2 ) + p (u 2 + ^ 2 + d d )+d pli 2 v + d oj 2 v ] 

c cr cr c c r rr cc rr c c cc r c 

+k [ - B (g 2 - d 2 ) - 6 (. 2 + .. 2 + d d )] 

c cr cr cc r rrcc 

2 2 

+ «.(-d xt - d • B ] 
rr c c cc r c 

2 2 2 2 2 2 2 
+ . (g -d ) + .(. + . +d d )+d d w 

c cr cr cc r rrcc ccrrc 


+ . + d ~ . 

c r c c r 


hk (, d r p + x v - d bp -.tv) 
cc r r err rr c c r c c 


+h•(. Bp - ■ : P ) 

err r c c 
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2 2 2 2 
+h(d uip + d o> p + u> u v ) 
rr c c cc r c c r c 


+k(-d w^B - d o^B ) - ) 

rr c c cc r c rcc 


4 2 2 

+d w + 2d u oj 

rr c cc r c 


M_ = hk(a) 6 p 
0 err 


p ) + hoj w p - kaj^cj^B + oj u) 

rcc ere rcc re 


Forming the Routh array for a general 6th-order polynominal, 


where 


(MM - M M ) /M 

5 4 6 3" 5 


(M 5 M 2 - m 6 m 1 )/m 5 


M(MM -MM)- 
3' 5 4 6 3 

M 5 (M 5 M, 


6 r 


<« 5 m 4 - m 6 m 3 ) 


M (M M -MM) 
V 5 4 6 3 ; 


H 5 H 0 


(MM -MM) 
k 5 4 6 3' 


P. = 


MM -MM 
5 2 6 1 


(M 1 (M 5 M 4 ' W " M 0 M 5 )(M 5 M 4 - W 
M 5 (M 3 (M 5 M 4 - M 6 M 3 ) - M 5 (M M 2 - MM)] 


( p i ° 2 - M 0 O 1 )/P 1 
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The most simple special case assumes 


1. No gyroscopic effects g = 0 

2. No off-diagonal damping, d = d =0 

cr rc 

3. Velocity sensing only p^ = p = 0 

4. Feedback only on velocity estimate, k = 0 
In this case 


M, = 1 


-6 l 4 2d + d + v h 

c cc rr c 


M = hSl (0 v - 8 v ) - £3 (d +d )+hv(d +d ) 


r r c c 


+d (d +d)+dd + 2u + 
cc cc rr cc rr c 


hi(d 6 v 
cc r r 


2 2 

d 3 v ) + hv (w + oj + d d ) 
rr c c c c r rr cc 


-2.8 (o> 2 + + d d ) + 2(d + d )J + d w + 2d u) 2 

c c r rr cc rr cc c rr c cc r 

2 2 2 2 
hi(u 6v - w 6 v ) + h(d w v + d w v ) 
err rcc rrcc ccrc 


-i(d <J0 2 0 + d a J 2 8 ) 

rr c c cc r c 


222 22222 
+w (uj +a.i + d d ) + d d ot + too) + d o) 

c c r rr cc cc rr c c r cc r 


. 2 2 2 2 4 „, 22 

hoi o) v - leu o p + d u + 2d to u 

crc rcc rrc ccrc 


2 4 

0) 10 

r c 


It is still impractical to work out the entire Routh array analytically, but 
we note that 


M. = l - 0 (first column must be >0 for stability) 
6 

Mq • 0 (last element of first column is positive) 

M r 0 implies v h - 8 ? > - (2d + d ) 

5 c c cc rr 
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N 1 M 5 


MM -MM 
5 4 6 3 


= h^JLv (3 v - 3 -v ) + h 2 v 2 (d + d ) 
err c c c cc rr 


-hi 3(3v - 3 v ) + h£(-43 vd - 23 v a +3vd 

err cc cccc ccrr rrcc 


+ 3 v d ) 
r r rr 


+hv [4d d + 3d 2 + d 2 + uj 2 ] 

c cc rr cc rr c 


+i 2 3 2 (d + d ) + (23 [-Ad d - 2d 2 - d 2 - w 2 ] 

c cc rr c cc rr cc rr c 


+ [4d 2 d + 2d d 2 + 2d 3 + 2d w 2 + 2d t/ + u) 2 d ] 

cc rr cc rr cc cc c cc r r rr 


This shows that the third entry in the first column of the Routh array will 
be positive under the more restrictive assumptions that 

h > 0 

l < 0 

B £ v - 3 r V f > 0 (more "control authority" on critical than 
residual modes) 

At this point, our algebraic capacity was exhausted, although it is clear that 
additional conditions can be derived. 

Conclusion: Since in practice d , d are very small, we have in ef- 

- cc rr 

feet shown that all but three conditions of Routh's test can be satisfied 
under conditions of a typical example, by choosing l < 0, h > 0. It is un¬ 
likely that the remaining conditions will allow any less-restrictive condi¬ 
tion than this. 
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SECTION 4 


AGGREGATE MODEL DEVELOPMENT 


4.0 Abstract 


A common practice in obtaining a reduced-order model of a large space 
structure for designing its vibration controllers is to truncate the 
vibration modes directly, according to, say, the line-of-sight motions. 

Such design models do not account for the influences of actuators on the 
truncated modes nor the influences of the truncated modes on the sensors. 

The attendant control spillover and observation spillover can cause the 
vibration controllers to destabilize the structures, even when these 
controllers have been carefully and satisfactorily designed to ensure the 
closed-loop stability with respect to the reduced-order design models. 

The development of aggregate design models was thus intended to provide 
better reduced-order models of the structures for use in the design of 
vibration controllers. Vibration modes can be aggregated with respect to 
actuator and sensor influences by a nonsingular transformation before 
making any truncation; the truncated portion can then be independent of the 
actuators and sensors. Such an aggregation transformation in the general 
form was carefully studied (see Subsection 4.1). A useful set of necessary 
and sufficient conditions and two separate sets of sufficient conditions 
for the aggregation transformation to be nonsingular were established; 
such conditions are useful as a guide in choosing specific values of a 
diagonal block of free parameters in the aggregation transformation. 


The aggregation transformation is not unique. Several interesting 
ideas of utilizing the uniqueness (i.e., of selecting the specific values 
of various free parameters in the transformation) for imposing additional 
properties on the aggregate design models were examined (see Subsection 4.3); 
some of their potential difficulties in obtaining closed-form analytical 
results were analyzed. 

Vibration modes can be aggregated also with respect to their 
influences on performance and the influences of disturbances in addition to 
actuator and sensor influences. For example, the control systems for a 
large space structure might be designed to minimize the transmission from 
the disturbance sources to the critical performance variables besides 
providing the structures with closed-loop stability. By including the 
influences of the disturbances as a "part" of actuator influences and, 
similarly, the influences on the performance as a "part" of sensor influences, 
this "generalized" aggregation problem was able to be recast into the same 
initial form (see Subsection 4.3), so that the same arguments and the same 
techniques could be applied and extended directly. 
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4.1 


Aggregation of Vibration Modes with Respect to Actuator and 
Sensor Influences 


A standard approach to model reduction is to truncate the open-loop 
modal representation directly, by neglecting (generally high-frequency) non- 
critical natural modes. This generally results in control energy being coupled 
into the neglected portion in a way which can be destabilizing. Here we explore 
a different choice of basis. Rather than using the open-loop modes directly, 
we aggregate them first based on the influence of the actuators and sensors. 

The advantage is that the coupling with the neglected portion becomes indepen¬ 
dent of feedback control gains, and is thus much more readily predicted. 
Furthermore, this choice of basis is also not affected by feedback. 

Consider the following usual finite-element modal representation of LSS. 

( n + ft 2 n = $ T b„u 

I F 

| y - c D tn + c v »n (4-1) 


where r) denotes an L-vector of normal coordinates; ft = {w^, .... ai^} and $ = 

[ d> , ..., 4> l ] are L x L matrices; and 4^ denote the natural frequency and 

the natural shape of mode i, respectively. L is usually a very large integer, 
n = n^) denotes the corresponding normal coordinates, u = (u^, .... 

u ) denotes the control inputs to m force actuators, and B„ the L * m influence 
m r 

matrix, y = (y^, ..., y^) denotes the observation outputs from displacement 
sensors and l^ velocity sensors (0 £ i, 0 <_ i); and denote the 

corresponding 1 x l influence matrices. Superscript T denotes transpose. Nor¬ 
mally, m < L and l < 2L. 

A corresponding state-space form is 



C = fC D <I> C v 4>] 
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We assume that both B and C are of full rank. The rows of matrix B and the 
columns of matrix C can be rearranged such that 


B = 


6 • [s i 


(4-3) 


with B^ being an n * m matrix of rank m and an l x n matrix of rank i for 

the smallest integer n. We assume the corresponding rearrangements have been 
made on the rows and columns of matrix A. 

Now, to aggregate the state variables of (4-2), namely, the vibration 
modes of (4-1), with respect to actuator and sensor influences, we aggregate 
the rows of matrix B and the columns of C by a similarity transformation T 
such that 


TB 


L 0 


CT _1 = [C,; 0] 


(4-4) 


where B^ is also an n * m matrix of rank m and also an l * n matrix of rank 
l. By defining 

x = Tx, A = TAT _1 , B - TB, C = CT _1 (4-5) 

we now have the following special but equivalent representation of the LSS 


| x =* Ax + Bu 
| y * Cx 


Partitioning these matrices conformably with (4-4) yields 


A = 


A 11 

A 12 

B = 

V 

— — 


A 21 

A 22 


_ 0 . 


lc 1 0] 


(4-6) 


where , in particular, is an n * n matrix. The corresponding state equation 
can thus be rewritten as 
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*1 = A 11 X 1 + A 12 X 2 + B 1 U 


X 2 A 21 x l + A 22 X 2 


(4-7) 


y = C.x 


11 


Using any such aggregation transformation (and^the orthogonal trans¬ 
formation representing the rearrangement of matrices A, B, and C before ag¬ 
gregation), we have aggregated the vibration modes with respect to the influ¬ 
ence of actuators and sensors. Observe, specifically, that the control of (4-7) 
does not directly affect nor is the observation directly affected by it. 

Neglecting we get the following aggregate design model of order only n: 


x = A^x + BjU 


y - C.x 


(4-8) 


Roughly speaking, such a reduced-order model represents an equivalent 
dynamics of those linear combinations of natural modes which can be directly 
measured and controlled. It is not an exact aggregation [1] - [3] of system 
(4-2) in the sense of Aoki. By proper selection of the aggregation transforma¬ 
tion T, this reduced-model may be made to approximate open-loop performance 
of system (4-1), much closer than any one obtained by direct truncation of the 
vibration modes. 


Comment 1 : In general, it is sufficient to have Z + m rows in the block B^ and 
S, + m columns in the block C^; hence n <_ £. + m. Depending on the specific 

form of matrices B and C at hand have, or due to the existence of both velocity 
and displacement measurements, the blocks B^ and need only max{£,m} rows 

and columns respectively; hence n £max{&,m}. Note that the possibility of 
having an aggregate design model (4-8) of order less than Z + m is an improve¬ 
ment over the preliminary results [7]. 


Comment 2 : The rearrangement of matrices B and C to have them bear the form 
described by (4-3) can be carried out using common computational procedures such 
as Gaussian elimination or Gram-Schmidt orthogonalization on rows of B and 
columns of C. Any prior effort in the placement of actuators (and sensors), 
or in the synthesis of actuator (and sensor) influences for preventing control 
spillover to (and observation spillover from) secondary modes, as discussed in 
[4-6], will facilitate the rearrangement of these matrices, since there will 
be fewer nonzero rows in B (and fewer nonzero columns in C) to work with. 

Since the number of actuator-influenced (and sensor-influenced) modes, partic¬ 
ularly those which are secondary modes, will be decreased, such spillover 
prevention efforts will also facilitate the aggregation and reduce the coupling 
Aj£ and A 21' 
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Comment 3 : The aggregation of, say, matrix B as described in (4-4) is equiv¬ 
alent to making proper linear combinations of the rows of B to produce the 
specified form. In conjunction with representation (4-7), this can be inter¬ 
preted as synthesizing the actuator influences so that control inputs will not 
spill over part X 2 of the aggregated state x. Compare this interpretation with 

the idea [4-6] of preventing control spillover to secondary modes by synthesis 
of actuator influences. To prevent control spillover by synthesis is to make 
proper linear combinations of columns of matrix B to achieve a partitioning 
similar to what is described by (4-4), namely, to determine a transformation T 
such that 


Br 



It follows from this comparison that the technique developed for computing 
synthesis transformation T can be modified to help compute the aggregation 
transformation T. This comment is equally applicable to the aggregation of 

matrix C with respect to an aggregation transformation S = T observation 
spillover, and synthesis of sensor influences. 

Let the transformation be partitioned with respect to (4-3) and (4-4) 
as follows 



The dimension of the blocks are as follows: 

T u , n x n; T^, n x (2L - n); T^, (2L - n) x n ; T^, (2L - n) x (2L - n). 
Then Eq. (4-4) can be written as 


T 11 B 1 + T 12 B 2 

" B 1 

(4-9) 

T 21 B 1 + T 22 B 2 

= 0 

(4-10) 

C 1 T 11 = 

6 1 

(4-11) 

C 1 T 12 = 


(4-12) 

is of the full row rank, its right 

inverse C B exists. 

Therefore, from 


(4-12), we have 


52 








( 4 - 13 ) 


T 12 ‘ C lS 


where 


R T T — 1 

C 1 - C 1 <C 1 C 1 ) +N 


12 


and is an arbitrary matrix such that 


C 1 N 12 = 0 


Similarly, has its left inverse and, from (4-10), we have 


T 21 T 22 B 2 B 1 


where 


“i E (b^i)-‘b^ + n 21 


and N 21 is an arbitrary matrix such that 


N 21 B 1 “ ° 


It follows from (4-14) and (4-16) that 


T = 


r t. . c R c„' 


'i o' 


11 12 




-T„„B„B^ T„„ 


0 T„„ 


22 2 1 22 


22 

L 


11 

* 2 B 1 


Therefore, the aggregation transformation T is nonsingular if 
T 22 and |t h + are nonsingular. 

Comment 4 : The specific value of matrices and B^ in (4-4) 

should be consequences of the aggregation transformation T. 
contrary that C. has been prescribed before the aggregation. 
(4-11) 


11 


c R c 

c ri 


(4-14a) 


(4-14b) 


(4-15) 


(4-16a) 


(4-16b) 


and only if both 

, in general. 

Suppose on the 
Then, from 
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Consequently, the matrix sum 

T U tC lWl ‘ C l( e i + e 2 6 2 S l) 

cannot be nonsingular unless for such a very special case that n = l >_ m and 
Cj + ^^B^ is nonsingular. Therefore, it is better not to prescribe the 
value of before aggregation. 

Matrix should not be prescribed, either. Premultiplying (4-9) by Cj^ 
and using (4-11) - (4-12) and (4-3), we have 

C 1 B 1 = + C 2 B 2 H CB (4-17) 


This relationship in general will not be satisfied if B is assigned a specific 
value before aggregation. 

In the sequel, we assume that both T^ and T 2 2 are nonsingular, and that 

both C^ and B^ are computed after aggregation. Now, from (4-13) and (4-15), we 
have 


T = 


11 

0 


22 


T~ 

x irr2 




Theorem 1 ; The aggregation matrix T is nonsingular if and only if the 
n x n matrix |l + T 1 |c^C 2 B 2 B^j or the (2L - n) x (2L - n) matrix 


' ~ _ T _1 D~ 

I + B 2 B^T^C^C 2 I is nonsingular, where 




‘ e i T n 


(4-18) 


Theorem 2 : The following are two separate sets of sufficient conditions 
for the matrix T to be nonsingular. 

(a) B 2 * 0 or C 2 = 0 

T ~T 

(b) T ii T n = <1 = ra, and C * EB for some nonsingular matrix E 
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Comment 5 : Condition (a) appears to be trivial, since either the actuator in¬ 
fluences or the sensor influences have been aggregated properly. This condi¬ 
tion, however, clearly indicates the advantage of using Gram-Schmidt or 
Gaussian elimination procedure for rearranging the rows of B and the columns 
of C, and the advantage of prior effort made in preventing control and obser¬ 
vation spillover by placement or synthesis; see Comment 2. 

Comment 6 : First part of condition (b) means that T is a special kind of non¬ 
singular matrix that is orthogonal; the identity matrix obviously is such an 
orthogonal matrix. The second and third parts mean that there are as many 
sensors as there are actuators and that the sensor influences as a whole are 
similar to the actuator influences as a whole. For LSS having only force 
actuators, this implies that no displacement sensors are considered for aggre¬ 
gation; namely, C p = 0 in (4-1). The special case in which E = kl represents 

the commonly considered case where each force actuator is collocated with a 
velocity sensor, and vice versa. Notice that vibration modes with such a 
special case can always be aggregated to achieve the desired representation (4-7). 

Comment 7 : The expressions (4-13) - (4-16) for T^ and T 2 i are slightly dif¬ 
ferent from those given in [7]. These slight improvements, however, have made 
it possible to obtain useful new results, such as Theorems 1 and 2. 

4.2 Ideas for Specifying the Nonunique Aggregation Transformation 

The similarity transformation T for accomplishing the aggregation as des¬ 
cribed by (4-4) is not unique; subject only to some rather minor constraints, 
submatrices T^, N j 2 > and N 2 1 are virtua lly free parameters. Additional 

properties may thus be imposed on the aggregate design model (4-8) by specifying 
these parameters. The following are some interesting ideas we made a quick 
look at. 


4.2.1 Minimization of Coupling Coefficients ,, and A ^^ 

Unlike the reduced-order design model usually obtained from direct 
truncation of vibration modes, the aggregate design model (4-8) has the desirable 
property that the coupling between the retained part x^ and the truncated part 

remains unchanged even in the closed loop with feedback control systems. 

This constancy does not depend on how the control systems are designed, be they 
designed as LQ regulators and observers, or as dynamic output feedback com¬ 
pensators [7]. Therefore, it would be very desirable if the open-loop coupling 
coefficients A^ and A£^ could be minimized by a suitable choice of the non¬ 
unique aggregation transformation. 

Partition the aggregated matrix A defined in (4-5) and the original 
matrix A conformably with the similarity transformation T. Then we have 
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r 


A 11 

A 12 

T 

11 

T 

12 


hi 

T 

12 

A 11 

A 12 _ 

1 - 

> 

ho 

A 22_ 

T 

21 

H 

ro 

l_ 


T 

.21 

T 

T 2 2 J 

A 21 

A 22_ 


Solving for and yields 

A 12 ' [ T U A 12 + T 12 A 22- (T 11 A 1, +T ,2*21> T n T 12] (T 22 

A 21 ‘ [ T 21 A 11 + T 22 A 21 * «21 A 12 + T 22 A 22 )T 22 T 2l]«11 


T 


T 


21 T 11 T 12 ) 

12 T 22 T 21 ) 


-1 


-1 


(4-19) 

(4-20) 


In terms of A^ and A^, we have 


11 


(T 11 A 11 + T 12 A 21 " A 12 T 21 )T 11 


A 22 = (T 21 A 22 + T 22 A 22 ~ A 21 T 12 )T 22 


Substituting (4-13) and (4-15) in (4-19) and (4-20), 


12 


21 


= T 


= T 


2 + T n c iS - A 11 T U C 1 6 2 


ll( A l 

T 11 C 1^2 A 21 T 11 C 1^2) ( X + ®2®1 T 11 C 1 C 2) T 22 

22 ( _ ® 2 ® 1^11 + A 21 “ B 2 B 1 A 12 B 2 B 1 
+ A 22 B 2 B i)( I + T n c ^2 B 2 B 2)~ lT n 


(4-21) 

(4-22) 


(4-23) 


(4-24) 


With the expressions (4-14a) and (4-18), the terms within the parentheses of 
(4-23) and (4-24) are virtually independent of and T^. Thus, to minimize 

A^ i- s to minimize and maximize T^> whereas it is the opposite for mini¬ 
mizing Therefore, the coupling coefficients A^ and A£j cannot be mini¬ 
mized simultaneously by any choice of T^ and ^ 22 ' 
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The only possibility lies in the selection of the null-space matrices 

£ 

and involved respectively in the right inverse and the left inverse 

B^. The terms in the first parenthesis of (4-23) can be minimized by a separate 
choice of N^« while those of (4-24) can be minimized simultaneously by a separ¬ 
ate choice of N,^. The matrix products in the second parenthesis of (4-23) and 
(4-24) can be maximized simultaneously by a joint choice of and 

The constancy of the terms in the parentheses of both (4-23) and (4-24) 
can be seen more easily if we consider the class of orthogonal matrices for 
T^, namely 


For this class 


T~V R r 


= C R C 
i z 


(4-25) 


which is independent of T^. All the terms within these parentheses are 
completely independent of T^! 

Nevertheless, the numbers of nonzero columns in and nonzero rows in 
A^^ can be reduced. Yeung Yam of HIT, while expanding the preliminary results 

in [7], has recently derived a specific choice of the free parameters in T 
such that, for the case n = £ + m, only n columns of A^ and n rows of A^j 
are nonzero. 

4.2.2 Minimization of Differences in Closed-Loop Transfer Functions 

Feedback control systems are to be designed based only on the reduced- 
order aggregate design model (4-8). The resultant closed-loop transfer functions 
for the design model alone generally will be different from the resultant ones 
for the overall system representing the LSS. It would be very desirable if 
the differences could be minimized by suitably specifying the free parameters 
of the aggregation transformation. 

As in [7], consider fixed-order dynamic compensators of the form 


| z = Fz + Gy 

| u = Hz + Ky + u (4-26) 

where the compensator state, z has a fixed dimension which is usually much 
less than 2L; u denotes the command inputs; F, G, H, and K are design 
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parameters to be determined based on the reduced-order model (4-8). Substitu¬ 
ting (4-26) in (4-8) yields the following closed-loop system of the design model. 


where 


2 


r z’ 



, 

l 

“ A n 

N 

+ 

M 



L J 



(4-27) 


11 


BjH 


GC, 


A 11 + B 1 KC 1 


(4-28) 


The transfer function matrix of the closed loop system (4-27) between u and y 
is found to be 

Q(s) = [0 C ](sl - A ) -1 ro' 

L b iJ 

- cjsl - A n - BjKCj - BjH(sI - F)' 1 GcJ " 1 Bj (4-29) 


Now substituting (4-26) in the overall system (4-7) yields the following 


closed-loop system. 


= A 


[0 C L 0] 


(4-30) 


where 


A = 

’ A ll_ 

_ A 12 


F 

B 1 H 

gc 1 ; 

A 11 + B l Kc li 

0 

A 12 


A 21 

A 22 


0 

A 21 ! 

A 22 - 
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with the same A^ as defined by (4-28). The tmaster function matrix of the 
closed loop system (4-30) between u and v is found to be 


1 _ 


Q(s) = C, si - A. - B KC - B H(si - K) c;C, - A 


1 1 


1 12 


, (si - A ) *A. I ' 1 


21 ; “1 


Define the difference by 


(4-32) 


E(s) = Q(s) - Q(s) 


(4-33) 


Then, one could attempt to choose the parameters T^, T 77 , N , anc ^ ^21 suc * 1 
that 


J 'trace E(s)E^(S)J dc. 


(4-34) 


is minimal, where s = o + jj and s denotes the complex conjugate of s. To ob¬ 
tain a solution to this minimization problem is impra;tical, if not impossible, 
judging from the algebraic complexity involved when the expressions (4-19) - 
(4-22) for A., and (4-13) - (4-16) for T.~ and T„. are substituted in (4-29) 
and (4-32). 1J 

Intuitively, the main difference between Q(s) and Q(s) is the term 
A^^I - A 22^ A 21” The minimization of this term will also minimize the over¬ 

all difference as defined by (4-33) or (4-34). See the preceding subsection 
(Subsection 4.2.1) for discussion on minimizing the coupling coefficients A^ 7 
and A 21> 

Another intuitive alternative is to maximize the terms within the brack¬ 
et of (4-29) so that the corresponding terms in (4-32) dominate the term 

A)2 (si “ A 22^ A 21‘ 

These two intuitive approaches may be much easier to take than direct 
minimization of the overall difference defined by (4-34) and any other way using 
(4-33). The necessary algebra still can be very complicated. 

4.2.3 Assurance of Closed-Loop Stability of LSS 

Again, feedback controllers are to be designed using only the aggregate 
design model (4-8), and at least to assure closed loop stability with respect 
to the design model. It would be desirable if closed loop stability could 
still be assured when such controllers were implemented on LSS, which contains 
not only the design model but also the truncated part. Such an assurance is 
generally impossible for the case of direct truncation of vibration modes. Re¬ 
call that the coupling A^ and A 21 of the aggregate design model are independent 

of any controllers designed by any method. Thus, if closed-loop stability of the 
LSS having a controller designed for the aggregate design model using some 
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model using some method can be guaranteed, then closed-loop stability of the 
LSS having any such controllers but designed by other methods can also be 
guaranteed. 


Consider again the feedback controllers of the form (4-26). Assume the 
design parameters F, G, H, and K are computed based on (4-8) so that the closed- 
loop system (4-27) is stable. The question is whether it is possible to find 
(and specify) among the nonunique aggregation transformation one such that the 
closed-loop system (4-30) is also stable. This requires the use of the stability 
theory for interconnected large scale subsystems; but the development of such 
stability theory has just begun [8] and no results are directly applicable to 
(or general enough for) our problem at hand yet. Moreover, with the expressions 

for A.., and T.. substituted in (4-31), the resultant matrix A is a very com- 
ij ij 

plicated nonlinear matrix function of parameters T^, T 22 ’ ^12* an< ^21' ® e_ 

riving an analytical stability condition for system (4-30) is very difficult, 
if not impossible, at present. An attempt was made, however, in applying 
Willems' generalization [9] of Lyapunov theory to system (4-30). 

4.3 Aggregation of Vibration Modes also with Respect to Tracking and 
Disturbance Rejection 

In the previous derivations, it was implicitly assumed that the pur¬ 
pose of the fixed-order compensator (4-26) was merely to enhance the stability 
of the closed-loop system. The aggregation was thus based only on the actu¬ 
ator and sensor influences. Now assume that there are a finite number of 
disturbance sources on the structure and a finite number of critical response 
variables (such as LOS motion in x and y directions). Then the compensator 
might be designed for minimizing the transmission from the disturbance sources 
to the critical response variables for all critical frequencies. This implies 
that the compensator places particular emphasis on damping out those vibration 
modes which can influence the critical response variables and can be excited 
by the disturbance sources. For such a generalized purpose, aggregation of 
the vibration modes should also include the influences of disturbance sources 
and critical response variables. 

We extend the preceding studies in the following way. Suppose the 
original representation is now expanded as follows: 

x = Ax + Bu + Dv 

y = Cx + Ew 

V r = Mx + Nu 


(4-35) 


bO 







where 


v is a vector of disturbance inputs 
w is a vector of observation noises 

r is a vector of response (not necessarily measurable) 
u, x, y are as before. 

Now define 


*0 

0 

E ' 


u 




» U = 

V 

_N 

0 

0_ 


- w _ 



r c i 


y 



1, y = 


C = 

i 




[m_ 


_ r _ 


x = &x + Bu 

y = Cx + Du (4-36) 


Now by arguing as before, there is a transformation which brings the repre- 

a — l - - -1 

sentation (4-35) into the form (A,B,C,D) = (TAT , TB, CT , D), where 



A il 

CM 


V 

A = 

A 21 

1 

CM 

CM 

C 

, B = 

0 


here Aj^ has the dimension between max {dim u, dim y} and (dim u + dim y). 

The representation (4-36) is now transformed into a similar form to (4-7). 
By making the same truncation, we obtain an aggregate design model which 
is also in a similar form to (4-8). 


Now suppose we have a compensator 


z = Fz + Gy 
u = Hz + Ky 
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Then the closed-loop dynamics is given by 


GC 


* * * 
B H A + B KC 


r = [NH M + NKC ] 


where N = N, E = E (since D = D), and B , C , D , M are identified by parti¬ 
tioning B and C, i.e.. 


2 


0 GE 


V 


+ 

* * 



X 


D B KEJ 

_ w - 





- 

z 



V 



+ 

[0 NKE] 



X 



-V 

- 

* 

* 

* * 




_ 


r * * 1 

B, 


B. D, 0 

1 

r * * i 

1 1 


in 

o 

o 


_0_ 


. 0 0 0. 


r- * r~ 'k 

C 


C = [c 1 0] = 


c 1 0 


0j 
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SECTION 5 


PARETO-OPTIMAL CONTROL 


5.0 Abstract 


Canavin's "modal-dashpot" concept and Aubrun's "low-authority 
control" theory were two off-the-shelf examples used to illustrate the 
multiple-objective optimization problems underlying the design of vibration 
controllers for large space structures (see Subsection 5.1). Canavin used 
pseudo-inverses to compute the output feedback gains required for augmenting 
a certain level of damping to primary modes (or, specifically, critical 
modes); the computed gains were of dismally high values. A careful exam¬ 
ination showed that such a high-gain problem need not exist because there 
were free parameters that could be adjusted to minimize the magnitude of the 
required gains. It also showed that the desire for large possible achievable 
damping and for small possible required gains could easily be formulated as 
a two-objective optimization problem and solved directly and efficiently by 
multiple-objective methods without computing the pseudo inverses and the 
appropriate values for the free parameters. In Aubrun's theory, the feed¬ 
back gains were to minimize a cost functional that was the sum of two 
disparate design objectives: the magnitude of the errors between predicted 
and desired closed-loop damping, and the magnitude of feedback gains. The 
tradeoff between the closed-loop damping errors and the feedback gains, and 
the variation of the former versus the latter, could be systematically and 
efficiently made by formulating the underlying two-objective optimization 
problem and then applying a multiple-objective optimization method. 

The applicability and application of two multiple-objective opti¬ 
mization methods were demonstrated by Pareto-optimal designs of modal 
dashpots for the 12-mode Draper model (see Subsection 5.2). Not only a 
sequence of the two-objective Pareto optimal designs but also the maximum 
damping achievable for each level of feedback gains were computed system¬ 
atically. Insights into the movement of the closed-loop poles as the 
feedback gains increase their magnitude were easy to obtain from the resulting 
root loci. Systematic evaluation for robustness of the resulting sequence 
of designs to parameter variations and to spillover became easy. The 
variations of the sensitivity (i.e., of the robustness) were also obtained. 
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5.1 Applicability of Multiple-Objective Optimization 

Methods to the Design of Vibration Controllers 

The design of a control system normally has more than one performance 
objective to achieve: e.g., large modal damping, low control response, little 
small feedback gains, small errors or overshoots, fast response, little 
spillover. As a standard practice, a designer usually combines the multiple 
objectives into one single performance index using some rather arbitrary 
weights. The designer, usually not satisfied with the resultant control 
systems, then repeatedly, aimlessly, and sometimes endlessly, changes the 
weighting coefficients to obtain new designs. Such frustration is not 
uncommon in the design of "modern control systems" using the linear- 
quadratic regulator theory; the weighting matrices Q and R frequently are 
changed over and over again. Experience with artificial single-objective 
optimization used on multiple-objective design problems has often left the 
designer dissatisfied with the lack of information (or systematic information) 
on the available tradeoffs between his multiple objectives and his "multiple- 
objective optimal" designs. Realistically, incompatible objectives cannot 
be unified; important multiple objectives should be addressed directly and 
their trade-off be made efficiently and systematically. 

Pareto optimality is a practical notion of "multiple-objective" 
optimality that is useful to control systems designers faced with multiple 
disparate performance objectives. Stated simply, a design is Pareto optimal 
if none of the multiple objectives can be further improved without degrading 
any other. Multiple-objective optimization by the method of proper equality 
constraints (PEC) or the method of proper inequality constraints (PIC) help 
the designers obtain both Pareto optimal designs of control systems and 
efficient tradeoffs of the multiple objectives directly and systematically. 

Consider the following two different approaches to the design of 
vibration controllers for large space structures: Canavin's "modal- 
dashpot" concept [1-4] and Aubrun's "low authority control" theory [5,6]. 

Let us point out their underlying multiple-objective design problems and 
briefly demonstrate the applicability of the above-mentioned multiple- 
objective methods [7-10]. 

Assume that large space structures are represented in modal coordinates 
as follows: 


n + 


n 2 n 


T 

= $ 


(5-1) 


y = c v <J>n 


(5-2) 
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where ft = diag joj^,... ,oj^ }, and $ * [4^,...,$^] are L*L matrices: and 

4>^ denote respectively the natural frequency and shape of the ith normal 

mode of vibration, n = (q^.q ) denotes the corresponding normal 

coordinates, u = (u,,...,u ) denotes the control inputs to m force actu- 
1 m 

ators, and B denotes the Lxm influence matrix, y = (yy.) denotes 
F IX- 

the observation outputs from l velocity sensors: denotes the correspond¬ 

ing £xL influence matrix. Superscript T denotes transpose. 


In general, a very large number of vibration modes are included in 
this model. Let {u) pi> * = 1»***>P» denote a selected subset for 

suppression; we call them primary modes for a self-explanatory reason. 
Completely neglecting all nonprimary modes yields the following truncated 
form 


fi p + ^n p 
y 


♦Jv 


Vp n p 


(5-3) 

(5-4) 


where ftp and <t p are similarly defined in terms of the p primary modes, 

and n p = (l pl ,,,.jhpp) denotes the corresponding normal coordinates. 

Vibration control systems are then designed for this model as if it had 
exactly modeled the LSS in question. 

Both Canavin's and Aubrun's approaches are for designing constant-gain 
velocity-output feedback control in the form 

u = -G v y (5-5) 


Substituting (5-5) in (5-3) and (5-4) yields the following closed-loop 
system dynamics 

n p + $ P B F G V C vVp + = o (5-6) 


Consider Canavin's modal-dashpot concept first. In order that each 
primary mode may be damped independently and arbitrarily, the feedback gain 
matrix G v is computed using pseudo-inverses as follows: 


■ K b f) i [*p b f(*X) t ]' 1 <2 w [(Vp) 1 Vp]' 1 (sM 


(5-7) 
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where Z p denotes a diagonal matrix of desired damping ratios £ for the p 

primary modes. It is assumed that S, > p and ra i p. Each actuator was 
assumed to be collocated with a velocity sensor in his numerical 
example [1,2]. For 10% damping for each primary mode, the computation using 
(5-7) resulted in a gain matrix of shockingly large values. He felt that 
"modal dashpots may be of limited utility due to the high gains produced 
by this approach" [1,2]. 

Such a high-gain problem actually need not exist. For exactly the 
same purpose of independently and arbitrarily damping the primary modes, 

T 

the gain matrix G can be computed using a right inverse of $ B^, and a left 
V r r 

inverse of C v <I> p , instead, as follows: 


°v - (*X) R ( 2z A)( c v" , p) L 


\R _ 

/ ,T \T 

f\T„ / 

' T \T "] 

) = 

(Vf) 

KM 

Vf) J 


( C vM " [( c vV) Vp ] ( c vV) + N v 


with N p and N v being arbitrary null-space matrices such that 


Wf = 0 


Wp = 0 


(5-8) 

(5-9) 

(5-10) 

(5-11) 

(5-12) 


(See Section 6.2, particularly Comments 7 and 8, concerning a generalization 
and an extension of the modal dashpot concept.) The similarity and the 
difference between (5-7) and (5-8) are obvious. Whenever the rank of 
T 

<I> B„ is less than the number m of actuators, matrix N_ need not be zero. 

P F f 

Similarly, whenever the rank of C y^p is less than the number l of sensors, 
matrix N„ need not be zero, either. Nonzero N and N serve as free para- 
meters and can be used to reduce the magnitude of gain matrix G^. The unde¬ 
sirable consequence of neglecting these adjustable free parameters is thus clear. 


Now that we have recognized the existence and advantage of free 
parameters N p and N^, an obvious question is how to choose their values 

to minimize the magnitude of the required feedback gains while using the 
feedback gains to augment as much damping to the primary modes as possible. 
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The underlying two performance objectives now become evident, and are 
formulated as a bi-objective optimization problem as follows. Note that 
in the interest of improving model fidelity, we also retain in the control 
design model some other important modes which we call secondary modes and 
denote by subscript S. Active damping of secondary modes is not required, 
but prevention of control spillover to and observation spillover from them 
is desired. Both primary and secondary modes are altogether referred to 
as modeled modes and denoted by subscript H. Therefore, the bi-objective 
modal-dashpot design problem is to find the feedback gains g^ and the 
damping ratios that 
maximize 

h ■ jb vpi (5 - 13) 

and minimize 

m SL 

J 2 = E E g* (5-14) 

i=l j=l 

subject to 


<f M B F G V C v' I> M 

2Z A 

(5-15) 

? i - ? Pmin 

for primary modes 

(5-16) 

? 1 - Smin 

for secondary modes 

(5-17) 


where weighting factors w^ are useful for expressing desired mode movement 

(see Subsection 5.2). Observe that no computation of the pseudo-inverses, 
the right or the left inverses, or the free parameters is required. They 
are all embedded in the equality constraint (5-16) of closed-loop modal 
decoupling. The method of proper equality constraints (PEC) is to set 
the first objective to a parametric level, i.e., 

£ w 5 = a (5-18) 

i=l 1 F1 


whereas the method of proper inequality constraints (PIC) is to bound the 
first objective by a parametric level, i.e., 


L w.e * a (5-19) 

i-1 1 Pi 
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Whether the method of PEC or PIC is used, the bi-objective problem (5-13) - 
(5-17) is all converted to a parametric single-objective quadratic 
programming problem of the following standard form: minimize 

J = x T Qx (5-20) 


subject to 


Ax ^ b 
x > 0 


Many algorithms and computer programs have been developed for solving the 
quadratic programming problem (5-20); see, for example, References [11-15]. 

For each selected value of the parameter a, the resulting quadratic 
programming problem is solved for the single-objective optimal solutions. 

Using the test for Pareto optimality as given in References [7-10], the 
Pareto optimal values of feedback gain g and damping ratios are 

obtained systematically in the increasing order of a. The resulting 
minimized values of the objective (i.e., the magnitude of gain matrix G^), 

plotted versus parameter a, will help the designer make optimal tradeoffs 
efficiently and systematically between achievable damping and minimum 
required gains. See Subsection 5.2 or Refs. [16-17] for more details. 

It is worth mentioning that the free parameters N and N are 

r V 

implicitly used to minimize the magnitude of required feedback gains and to 
d ecouple secondary modes from primary modes in the closed loop . Such was 
not possible with Canavin's origin modal-dasnpot concept. On the other hand, 
there is no question at all about the applicability of multiple-objective 
optimization methods to design of modal dashpots. 

Now, consider Aubrun’s low-authority-control (LAC) theory. According 
to Aubrun [5,6], for sufficiently "small" feedback gains g„ , the closed- 

loop damping of mode i is predicted using the following formula [5,6]: 

dA. * 2C..U). * ^B F G vVi (5 " 21) 

using the open-loop frequency oo^ and mode shape <Jk. To compute the gains g^, 
let (dA^)p denote the predicted closed-loop damping given by the above 
formula, and let (dA^p denote the desired closed-loop damping imposed by 
the LAC controller design. Then the gains g„ are chosen so as to minimize 
the quadratic cost functional 

J(C V ) - E» t [ «!Vd - (d Vp ] 2 + £ 8 ij <5 - 22) 
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in which the modal weights w^ help specify pole locations and the 

term improve robustness of the controller [6]. It is not difficult 

to see from (5-22) that two disparate objectives, one concerning damping 
and the other concerning gains, are "unified". The underlying bi-objective 
optimization problem is to 
minimize 


J 1 * (d Vp] 2 


(5-23) 


and minimize 


J 


2 



(5-23) 


Pareto optimal solutions of this problem by the method of PEC or PIC will 
provide a curve of J- versus for systematically making efficient tradeoff 

between achievable damping errors and the corresponding feedback gains 
required. The first objective, J^, is converted to a parametric equality 

constraint if the method of PEC is used, while it is converted to a 
parametric inequality constraint if the method of PIC is used instead. All 
converts the bi-objective problem to a simplest form of nonlinear programming 
problem, where the single objective function and the parametric constraint 
are quadratic. Many nonlinear programming algorithms and computer programs 
already developed [14,15] can be used in solving the resultant constrained 
optimization problem. Applicability of the multiple-objective optimization 
methods is, again, not a problem at all. 

Comment 1 : Both underlying bi-objective design problems (5-13) - 
(5-17) and (5-23) - (5-24), have the minimization of feedback gains as an 
objective; such is a very practical engineering objective in addition to 
the design objective on damping. 


Comment 2 : The closed-loop damping expressed by formula (5-21) may 
not be valid since the closed-loop modes may not remain decoupled with the 
same open-loop frequencies, even with respect to the reduced design model. 
Only when the gains g^ are extremely small will the closed-loop modes be 

virtually uncoupled. On the other hand, the desired independence can be 
explicitly imposed, such as (5-15); the resultant Pareto optimal designs 
then will certainly ensure that the open-loop modes remained decoupled at 
least in the closed loop of the design model. 


Comment 3 : The minimization of feedback gains, such as expressed by 
(5-14) or (5-24), does not actually provide robustness of the controllers 
designed. It will reduce control and observation spillover, however, 
because the residual modes will be less excited and less observed. 
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Application to Draper Tetrahedral Model of LSS 


The 12-mode Draper tetrahedral model of LSS was used as an 
example [16,17]. Pareto optimal designs of modal dashpots and the tradeoff 
curve were sequentially generated by a bi-objective optimization program 
developed. In search of more or less uniform movement of the primary modes 
under efficient tradeoffs of damping and gain objectives, we found the 
form of inverse-sum-of-squared-influence (ISS1) for the weighting factors w^ 
in (5-13) more satisfactory than others. 

Each of the resulting designs was then examined for its robustness 
to model parameter variation and spillover (resulting from modal truncation). 
For all designs with velocity sensors collocated with force actuators, 
robustness of stability to parameter variation and to spillover was clearly 
demonstrated. The damping performance is relatively insensitive to 
parameter variation and to spillover; however, the sensitivity rapidly 
increases as the desired level of critical-mode damping increases. In other 
words, modal dashpots can theoretically be designed to provide any desired 
level of damping to the critical vibration modes. In practice, however, 
only low designed levels of modal damping can be closely achieved due to 
this sensitivity to parameter variation and spillover. 

For those designs where the sensors and actuators were not all 
collocated, stability and damping performance are also relatively insen¬ 
sitive to parameter variation. However, spillover causes all designs to 
destabilize the originally stable structure; at least half of the resultant 
closed-loop poles have positive real parts. Destabilization by modal 
dashpots, which are also referred to as robust controllers, has not been 
demonstrated before. Consequently, collocation of sensors with actuators 
is essential to the robustness of modal dashpots when used as vibration 
dampers for realistic large space structures. 
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SECTION 6 


GENERAL REQUIREMENTS FOR STABILIZATION 


6.0 Abstract 

Feedback from velocity sensor outputs can augment damping and 
feedback from displacement sensor outputs can augment stiffness to large 
space structures. The design of the feedback gains are normally based only 
on modeled modes; the existence of unmodeled modes is completely ignored. 

Due to the attendant control and observation spillover, even stability of 
the closed-loop system having feedback controllers thus designed may not be 
ensured, let alone •'he additional damping and stiffness desired. Useful 
general conditions in various useful forms were recently established (see 
Subsection 6.1) for guaranteeing closed-loop asymptotic stability. Such 
conditions are useful as a guide to the design of closed-loop stabilizing 
feedback gains. carge space structures considered were not limited to 
those having only stable elastic modes; rigid-body modes and unstable modes 
were not excluded. For those special structures having only stable elastic 
modes and rigid-body modes, useful simple conditions were also established. 
The conditions for robustness to parameter variation and to modal truncation 
were also established. 

Canavin's modal-dashpot concept was generalized and extended (see 
Subsection 6.2). The similar concept of "modal springs" was therefore 
proposed. The general formulas for computing the feedback gains of 
generalized modal dashpots and modal springs were derived. These formulas 
contain various adjustable free parameters. The closed-loop asymptotic 
stability ensurance and the robustness to parameter variation and to modal 
truncation of modal dashpots and modal springs were established, thus at 
least providing a theoretic foundation for Canavin's original concept. 
Although no analytic proofs for the closed-loop stability ensurance and 
robustness of modal dashpots were given by Cnnavin when he proDosed the 
intuitive concept, a numerical demonstration was made with a typical 
example of large space structures. 

Control spillover to and observation spillover from secondary modes 
can be prevented before the beginning of the controller design process by- 
proper placement of actuators and sensors or by proper synthesis of 
actuator and sensor influences. A general-purpose computer program is being 
developed for synthesizing actuator influences so that control spillover 
to secondary modes can be prevented and the desired control influences on 
primary modes can be obtained simultaneously. A recent theoretical analysis 
has led to the conclusion that such a computer program can be used to 
synthesize sensor influences by duality and even to determine the proper 
location of actuators and sensors (see Subsection 6.3). 
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6 . 1 General Conditions for Closed-Loop Stab 11iy .at: ion by 
Output Feedback Control 

Consider the active control of a generic large space structure (LSS) 
which may contain rigid-body modes, unstable and stable elastic modes. As 
usual, let the modal equations generated by a finite-element computer program 
(such as NASTRAN) be 


+ 


.T 

v 


f 


( 6 - 1 ) 


where n denotes the L (normal) ’modal coordinates and 1 the natural mode shape.-.; 
f denotes the m generalized forces; A and Z are diagonal matrices denoting, 
respectively, the modal damping and Lhe modal stiffness coefficients. In gen¬ 
eral, the diagonal elements of A and Z may be positive, zero, or even negative, 
and can be extremely small if not zero. Note that if q denotes the generalized 
physical coordinates, then 


q = ; : (6-2 ,1 

Consider specifically the case of constant-gain output feedback control, namely, 
i velocity measurements and displacement measurements are fed back to m 

force actuators with constant gains as follows 


f 

= b f u 

(6-3) 

u 

= -"G y ~ G V T . 

(6-4) 


\rV DM) 


>’v 

= c v q 

(6-5a) 

y D 


(6-5b) 


In other words 


f 


B F°D C D q 


( 6 - 6 ) 


Substituting (6-6) in (6-1) yields 
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(6-7) 


where 


d a j b f g v c v 


k a 5 b f g d c d 


(6-8a) 


(6-8b) 


denote the additional damping and stiffness matrices respectively. The design 
of gain matrices and should make both the augmented damping matrix 

/ T \ / T \ 

l A + $ D <t> I and the augmented stiffness matrix l E + ^ symmetric and posi¬ 

tive definite, in addition to achieving certain design objectives (e.g., in¬ 
creasing the damping of critical modes to desired levels). The reason is sim¬ 
ple, since it is well known that a closed-loop system as described by (6-7) is 

asymptotically stable if both matrices /a + | and j are sym¬ 

metric and positive definite, and closed-loop stability is a basic requirement. 

Note, however, that the matrices A and Z are not positive definite un¬ 
less all their diagonal elements are positive. Also note that both matrices 
and are singular, since there are more modes represented by (6-1) than 

there are actuators or sensors placed on the structure (i.e., L •> i , L -■ ^). 
Moreover, the design of gains G^ and G^, in general, can be based only on a 

reduced-order model of the structure; in other words, onlv a truncated version 
of (6-1) can be used. Therefore, the asymptotic stability of the closed-loop 
system (6-7) does not follow automatically with an arbitrary design of G^ and G^. 

Let subscript M denote the part of (6-1) that is included in the re¬ 
duced-order model for control system design and subscript U the unmodeled 
(i.e., truncated) part. Then Eqs. (6-1) - (6-5) can be rewritten as follows 


1 + L M n M 

= 

T 

Vf 11 

(6-9a) 

I + V'u 

= 

T 

Vf u 

(6-9b) 

~Vv “ 1 

"D y D 

(6-10) 

C V <I 'M n M + 

V 

u n u 

(6-1 la) 

c pV'm + 

V 

u n u 

(6-1 lb) 


The design of gains G^ and G^ ignores the existence of the unmodeled part com¬ 
pletely, and is to make only the modeled part achieve desired damping and 
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stiffness objectives. Due to control spillover to and observation spillover 
from the unmodeled part, as can be seen from Eq. (6-9b) and (6-11), the overall 
closed-loop system (6-7) may not achieve the damping and stiffness objectives. 
In particular, the symmetry and positive definiteness of the overall matrices 

| A + ■?> ) and 1 j do not follow automatically from only making their 

modeled parts, j and j, symmetric and positive defi¬ 

nite. We have recently obtained the following general conditions (Theorems 1 
through 5) for guaranteeing the symmetry and positive definiteness of such 
overall matrices and hence for guaranteeing the asymptotic stability of closed- 
loop system (6-7). Such conditions are also useful as a guide in the design of 
gain matrices Gy and C Q . A useful special case is discussed in Subsection 6.2. 

Theorem 1 : Assume that is symmetric and | is positive 

definite. Then, for (a + 4>^B l ) t° be positive definite, it is necessary 
and sufficient that ' 


= A, 


IT A U 


T 

- <1 D <j 

U A 


M \ M 


..T_ t 

v m°a’m 


)"' 4 Svu 


( 6 - 12 ) 


is positive definite. 

Theorem 2 : Assume that K- A is symmetric and j is Positive 

definite. Then, for |l + to be positive definite, it necessary and 

and sufficient that 


'uVu ' 4 uYm\ ; m 


'M A'M 


-1 


T 

$ K i 

M A U 


(6-13) 


is positive definite. 

Comment 1 : For closed-loop asymptotic stability of the overall system, it is 
not enough just to guarantee the positive definiteness of the modeled parts 

( A M + Wm) anC * ( £ 'M + ^M^a'^m)’ According to Theorems 1 and 2, the matrices 
Aj and ^ defined by (6-12) and (6-1j), respectively must be made positive 
definite as well in the design of Gy and . This explains why current ap¬ 
proaches to the design of LSS controllers (e.g., direct applications of state- 
of-the-art control techniques only to the modeled part without any considera¬ 
tion of the unmodeled part) cannot guarantee closed-loop asymptotic stabilitv. 

Comment 2 : It is convenient as well as intuitively desirable to be able to 
find a design of, say, G such that the additional damping matrix for the 

V T 

modeled part is diagonal. For example, using the pseudo-inverses of and 

Cv'+m with a diagonal matrix will make I ^ diagonal. Notice, however that 

T 

those Gy designed using such pseudo-inverses will make only s y mmetr i c > to 




1 
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T 

make the whole matrix 4> D.l> symmetric requires that each velocity sensor be 
A T 

collocated with a force actuator such that C„ = k,,B for some nonzero constant 

V V F 

ky. Note, however, that collocation is not assumed in either Theorem 1 or 2. 

Comment 3 : These theorems are applicable to the general case of large space 
structures that contain not only stable elastic modes but also rigid-body 
modes and unstable modes. To provide a better insight into the requirements 
in stabilization of general large space structures, refine the partitioning 
of matrices A and 2 as follows. Denote the group by subscript G if their 
diagonal elements of A are greater than zero, and by double subscript LE, if 
otherwise. Similarly, denote the group'by G' if their diagonal elements of 
2 are greater than zero, and by L'E', if otherwise. Then we get the following 
similar general conditions 

Theorem 3: Assume that D, is symmetric and ( A,„ + is positive 

A \ LE LE A LL J 

definite. Then, for IA + to be positive definite, it is necessary and 

sufficient that 


T 

t> D <t 

G A 


; t d ,*> I 
G A LEI 


LE 


<p D 4> 

LE A LE 


r 


T 

0 D { 

LE A G 


(6-14) 


is positive definite. 

Theorem 4: Assume that K, is symmetric and ( 2, „ ,K. lT .,) is 

A \LE LLALE/ 

positive definite. Then, for I I + ) to be positive definite, it is 

necessary and sufficient that ' / 


'C' 


*G' K aV 


- $T g'Vl 


' E' (‘ 


L'E' 


+ t> 


L'E 


' k a*l' 


E' 


T 

♦l'e'V 


g' 

(6-15) 


is positive definite. 

Comment 4 : The unstable or undamped modes, having either nonpositive modal 
damping or nonpositive modal stiffness, naturally should be included in the 
model for designing gain matrices G^ 7 and G^, and at least both 

LE + "leVvVle) and ("L'E' + Ve' B F Wl’E' ) Sh ° uld be made P ositive 
definite. Theorems 3 and 4, however, show that this is not sufficient unless 
both an d ^ 2 ’ as defined by (6-14) and (6-15), are guaranteed to be positive 

definite. In general, flexible structures having small positive modal damping 
(i.e., small diagonal elements of A^), and lower natural vibration frequencies 

(i.e., small diagonal elements of 4^,) will be difficult in providing such 
guarantees. 

Consider the following common special cases: only rigid-body modes, and 
elastic modes, undamped and lightly damped, are present. In the above partition¬ 
ing, the groups denoted by LE and L'E' now correspond to only zero. 
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diagonal elements of A and £, respectively. Denoting these groups now by E 
and E', we get the following two very useful consequences of Theorems 3 and A, 
concerning closed-loop stability and its robustness. 

Theorem 5 : Assume that each unmodeled elastic mode has some positive 
damping, no matter how small it is. Then the closed-loop overall system (6-7) 

is asymptotically stable, i.e., both (A + {^D tland (z + <f T K *) are symmetric 
and positive definite if ' / V A / 


are symmetric 


(a) D. and K are symmetric and at least positive semidefinite, and 


(b) ^e^A^E an< ^ '* , E'^A' 1 E , are P os: *- t:: ‘- ve definite. 

Comment 5 : This theorem implies that the closed-loop asymptotic stability is 
robust to modal truncation, so long as conditions (a) and (b) are satisfied 
and all the unmodeled (i.e., the truncated) modes have some positive damping, 
finite or infinitesimal. First, unmodeled modes do not appear in conditions 
(a) and (b) at all, control spillover and observation spillover thus have no 
effect on closed-loop asymptotic stability so long as the assumption and con¬ 
ditions of Theorem 5 satisfied. Second, closed-loop asymptotic stability does 
not depend on how many damped elastic modes are included in the design model 
nor on which of them are, since neither conditions (a) nor (b) depend on any 
damped elastic mode. 

Theorem 6 : Suppose the gain matrices C; and Gy have been designed to 

satisfy conditions (a) and (b) of Theorem 5. Let there be parameter variations 
(or errors) in the mass matrix, the stiffness matrix, the natural frequencies, 
or the mode shapes of the structure concerned. Assume that 

(i) the number of rigid-body modes and that of undamped modes 
remains unchanged, and 


Where the overbars denote the presence of parameter variations (errors); C„ and 
C , are nonsingular matrices. Then the closed-lo-.p overall system remains 

Ej 

asymptotically stable. 

Comment 6 : Theorem 6 implies that closed-loop asymptotic stability is robust 
to changes in the shape of any damped elastic modes, modeled or not. 

The following subsection (Subsection 6.2) continues the discussion of 
this special case with a specific design of gain matrices and Gy that can 

provide not only the closed-loop asymptotic stability but also its robustness 
to modal truncation and parameter variation. Before concluding this subsection, 
we remark that the controllability and observability of all the rigid-body 
modes, all the unstable modes (with either negative modal damping or negative 
modal stiffness coefficients), and all the critical elastic modes, are assumed. 
The actuators and sensors should therefore be properly placed or their existent 
influences be properly synthesized, according to the techniques and conditions 
previously derived; see Refs. 1-3; Subsections 6.3 and 6.A. 
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6.2 Closed-Loop Stabilization with "Modal Dashpots" and "Modal Springs" 

A special output feedback control design technique is presented in this 
section for integrated control of elastic and rigid-body modes, for increasing 
both the damping and frequency of selected modes, and for providing closed- 
loop asymptotic stability and its robustness to modal truncation and to param¬ 
eter variation. This technique represents an improvement over previous work 
[4-8] concerning output feedback control of LSS or control of rigid-body modes, 
and is an advancement along Canavin's modal-dashpot design philosophy [4], 

Consider again the common special case of large space structures: 
rigid-body modes, undamped and lightly damped elastic modes may all be present. 
Assume that the rigid-body modes are not to be ignored as was in [4] but to 
be controlled along with the undamped elastic modes and those lightly damped 
elastic modes which are critical to the damping performance of the LSS. These 
modes constitute the primary modes. Assume that the design model is formed by 
these primary modes. Therefore, all other modes are unmodeled modes and all 
are assumed to have some positive damping, despite how small they are. Using 
the same partitioning as described at the end of the foregoing subsection, we 
have the following direct consequences of Theorem 5. 


Theorem 7 : Denote the number of modeled (i.e., the primary) modes by 
p. Assume that th^re are at least as many force actuators as there are modeled 
modes (i.e., m ^ p) , and that each force actuator is collocated with a velocity 

T T 

sensor and a displacement sensor, (i.e., Cy = kyBp, for some non¬ 

zero constants ky and k^). Assume further that all modeled modes are com¬ 
pletely controllable. Then the closed-loop system is asymptotically stable if 
the gain matrices Gy and G Q are computed as follows 



(6-16) 

(6-17) 

(6-18a) 

(6-19a) 

(6-20a) 


where A and £ denote two arbitrary positive definite p x p matrices desired 
M M 

as the additional damping and stiffness for the modeled modes; N , N , and N Q 
are any null-space matrices such that 
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$ m b f n f 

= 0 

(6-18b) 


Vv - 

4 

(6-19b) 


- 


(6-20b) 

Comment 7: 

It is technically convenient. 

though not necessary. 

to consider 


A„ and as diagonal matrices. Canavin [4] originally called the resulting 

modally decoupled velocity output feedback controllers "modal dashpots", in 
which the rigid-body modes were ignored, no displacement feedback was consider¬ 


ed, both Np and were zero, and A^ was diagonal. Here, we extend the concept 

in two ways. First, any velocity feedback controllers designed according to 
(6-16), (6-18) and (6-19) are called (generalized) modal dashpots, so long as 

is diagonal. Second, any displacement feedback controllers designed accord¬ 
ing to (6-17), (6-18) and (6-20) are similarly called "modal springs", so long as 

•k 

the matrix is diagonal. It should be noted that these terms (including the 

one initiated by Canavin) and, particularly, the implied modal independence in 
the closed loop make sense only with respect to the design model. In the 
presence of unmodeled modes, the modeled modes may not maintain their designed 
independence due to control and observation spillover. 

Comment 8 : Since the gain matrices <and given by (6-16) - (6-20) represent 

a special form that satisfies conditions (a) and (b) of Theorem 5, the asser¬ 
tions on robustness of closed-loop asymptotic stability to modal truncation 
and to parameter variation remain valid; see Comments 5 and 6, and Theorem 6. 

Comment 9 : Stability and robustness results for the special subcase [9] in 
which there are no undamped elastic modes, matrices N,,, N^, and in the design 

* ★ 

equations (6-18)-(6-20) are all zero, matrices A^ and Z^ in design equations 
(6-16)-(6-17) are all diagonal, and matrix C in Theorem 6 is diagonal were cb- 

tained jointly with Dr. Y.H. Lin of the Jet Propulsion Laboratory, California 
Institute of Technology (who was formerly with the C.S. Draper Laboratory). 
Various numerical demonstrations were also made, using a 10-meter free-free 
beam. 
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6.3 


Prevention of Control and Observation Spillover by Synthesis of 
Actuator and Sensor Influences 


Theorem 5 implies that for LSS having no unstable modes whatsoever and 
no undamped unmodeled modes, control and observation spillover will not cause 
closed-loop instability so far as output feedback control is concerned. The 
design of the output feedback controllers can be based entirely on the modeled 
modes which include all rigid-body modes, all undamped elastic modes, as well 
as critical elastic modes. Nevertheless, the closed-loop performance of the 
vibration controllers (e.g., the damping and frequency shifting of primary 
modes) will be degraded if control spillover and observation spillover exist. 
Theorems 4 and 5 imply that control and observation spillover problems may 
make it difficult to guarantee closed-loop stability of LSS having some unsta- 

T T 

ble modes, unless the control spillover matrices 1> B_ and 5 ,B , observation- 

U r ' L 


G' F 


spillover matrices 0^$^ and are all of small magnitude. It is now a 

standard understanding that control spillover and observation spillover will 
cause closed-loop instability for LSS using state-feedback controllers with 
state estimators. 


Control spillover to and observation spillover from secondary modes can 
be prevented by proper placement of actuators and sensors on the structures 
and by proper synthesis of actuator and sensor influences [10-12]. We have re¬ 
cently concluded that the technique proposed in [10-12] for synthesizing actu¬ 
ator influences can be extended not only to the synthesis of sensor influences 
(because of duality) but also to the placement of actuators and sensors. The 
formulas derived [10-12] for proper placement of actuators (and sensors) can 
be recasted into similar mathematical forms for proper synthesis of actuator 
(sensor) influences. A general-purpose computer program is being developed for 
synthesizing actuator influences so that control spillover to secondary modes 
is prevented and, moreover, desired control influences on primary modes is ob¬ 
tained . 


For preventing control spillover to secondary modes by synthesis of 

T 

actuator influences, the control spillover matrix LgBp is to be transformed 
using column operations Q and row operations P such that 


W 

w 


11 

W L2 

21 

W 22 


P(l'gB F )Q 


( 6 - 21 ) 


T 

with being nonsingular and having the same rank at 'l’gBp. Such transforma¬ 

tion (using Q and P) can be accomplished by Gaussian elimination, Gram—Schmidt 
orthogonalization, or singular-value decomposition. Preferably, W^ is made 

to be zero. Then the synthesizer T for prevention of control spillover to 
secondary modes, i.e., as a solution of 






0 


( 6 - 22 ) 


T 

W 


is given by 


' (Vn“u + « 2 ) f 


(6-23) 


where Q = [Qj» Q,, ]. The matrix F represents free parameters and can be adjust¬ 
ed for achieving^desired control influences on primary modes. Let B denote 
a desired control-influence matrix for the primary modes. Then F is to be 
solved from the following equations 

( _V 1 W 1 1 W 12 + V 2 )^ = ® (6-24) 

where 

[ V 1 ,V 2 j ■ t6 - 25 > 

The program under development is first to carry out the transformation 
(6-21) using Gaussian elimination, so that W = 0. The transformation (6-25) 

T * 2 

of actuator influence matrix 'l-B,. on primary modes is carried out simultan- 

r r 

eously, since the same operations Q are used. This program is then to compute 
f from (6-24). A desired synthesizer F is thus obtained from (6-23). Notice 
that Wj 2 in both (6-23) and (6-24) is zero. A very simple and specific version 
was developed for application to the Draper tetrahedral model of LSS. Modifica¬ 
tion, generalization, and addition are underway. 

For preventing control spillover to secondary modes by placement of 

actuators the corresponding formulas for determining the proper actuator- 

influence matrix B„ are as lollows 
F 


W,,W 0 ' = lip (6-26) 

I Z j o 

T 

with being nonsingular and having the rank of <J>g, which is equal to the 
number of secondary modes. B is to be a solution of 

r 


T 

<t> B_ 
S F 


= 0 


(6-27) 
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and is given by 


B F • K“iS + «2) B 


( 6 - 28 ) 


where Q = [Q^,Q^] and B represents free parameters. Let B denote a desired 
control-influence matrix for the primary modes. Then B is to be solved from 

(~ v l v l~ lv 2 + V 2) B = B (6-29) 

where 

[V v 2j - (6 - 30 > 

The one-to-one mathematical similarity between (6-26) - (6-30) and (6-21) - 
(6-25) is obvious. 

To place actuators to prevent control spillover to all nonprimgrv modes , 
if possible, is to find a matrix B such that the composite matrix Mlp§ is re- 

alizable as actuator influences, where M deontes the mass (inertia) matrix and 
the orthogonality of mode shapes is used. The influences of a real actuator on 
LSS are zero at some modes, and nonzero at others. So, the realization of 
the above ideal actuator influences by placement is to find B such that "con¬ 
trol spillover" to some modes is prevented while desired control influences at 
some other modes are obtained. The applicability and usefulness of the above- 
mentioned general-purpose computer program initially only for synthesis of 
actuator influences is obvious. 
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SECTION 7 


STOCHASTIC OUTPUT FEEDBACK CONTROL 


7.1 Introduction 


The minimum variance fixed form output feedback compensator has been 
recently proposed [1,2] as a possible solution to the problem of designing 
reduced order controllers for plants of large dimension. The zero order 
limit of this compensator, the stochastic output feedback controller [3], 
is of particular interest because it is the easiest such compensator to 
implement and simultaneously results in the simplest design problem. The 
price of this simplicity is that this compensator generally achieves poorer 
performance than optimal compensators of higher order. As a result, the 
stochastic output feedback controller should prove useful in assessing 
complexity/performance tradeoffs in the design of simple compensators for 
complex plants. 

Like the reduced order compensators considered earlier by Levine, 

Athans and Johnson [4,5,6], however, even in the zero order limit the 
minimum variance fixed form output feedback compensator gains and cost are 
difficult to compute, especially for plants of large dimension. Consequently, 
upper and lower bounds on the optimal constant gain cost that are relatively 
easy to compute have been developed and are described in the following 
sections. Specifically, Section 7.2 presents a more detailed discussion of 
the stochastic output feedback controller and its relation to several other 
compensator formulations such as the minimum order dual observer based 
compensator [1,7] and the Levine-Athans-Johnson reduced order compensator 
[4,5,6]. Section 7.3 derives an easily computed lower bound on the optimal 
stochastic output feedback controller cost while Section 7.4 derives three 
easily computed upper bounds on the cost. Finally, Section 7.5 discusses the 
principal results of the previous sections. 

7.2 Minimum Various Fixed Form Output Feedback Compensator 

Consider a plant of the form 

jc(t) = Ax(t) + Bju(t) + v(t) (7-1) 

with observations 

Z (t) = cx(t) + w(t) (7-2) 

Here x(t) £ R n is the plant state vector, jj(t) £ R m is the control input 
vector, y(t) e R^ is the observation vector and A, B, C are constant matrices 
of dimension nxn, nxm and p x n, respectively. The vectors v^(t) £ R n and 
w(t) £ are white noise processes representing plant driving noise and 





observation noise, respectively. These are independent stationary zero-mean 
processes with covariances 


and 


E{v(t) v'(T)} 
E{w(t) w'(T)} 


= V<S(t - x) 
= W6(t - t) 


(7-3) 


As defined in [1,2], the minimum variance fixed form output feedback compen¬ 
sator is of the form 


u(t) = Hz(t) + Mjr(t) (7-4) 

where _z(t) £ is the compensator state vector whose dynamics are described 
by 


£(t) = Fz(t) + Gjr(t) (7-5) 

Here F, G, h and M are constant matrices of dimension qxq, qxp, mxq and mxp, 
respectively. If mean-square stabilizing compensators of this form exist 
for given dynamic order q, the optimal compensator is defined as the one 
that minimizes the asymptotic state covariance weighted by the n*n matrix Q 
of design parameters with respect to the matrices F, G, H, and M, i.e., 

J = lira E{x'(t) Q x(t)} (7-6) 

r~K» 


Necessary conditions for the optimality of this compensator are presented 
in [1,2] for compensators of arbitrary dynamic order q. Two cases of 
particular interest are the stochastic output feedback compensator (SOFC), 
corresponding to q = 0, and the minimum order dual observer based compen¬ 
sator (M0D0BC), corresponding to q = n - m. 

For the SOFC problem, the necessary conditions reduce to 

(A + BMC)X + X(A + BMC)' + V + BMWK'B' = 0 (7-7) 

M = -(B’PB) _1 B’PXC'W _1 (7-8) 

P (A + BMC) + (A + BMC)' P + Q = 0 (7-9) 

assuming the required inverse exists. The n*n matrix X represents the 
asymptotic state covariance matrix so the cost may be expressed as 


J = tr[QX] 


(7-10) 
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As In the case of the compensators considered In [4,5,6] It is to be expected 
that the cost functional will exhibit local minimum [8] so th c equations 
(7-7), (7-8) and (7-9) will not possess a unique solution. 

Specific necessary conditions can also be developed for the MODOBC 
problem [1], Here the process is more difficult, but the solution can be 
shown to be unique [7] implying the conditions are also sufficient for 
optimality. These equations are derived by first transforming the problem 
into standard controllable form and then applying the necessary conditions 
for the minimum variance fixed form output feedback compensator. Determination 
of the optimal MODOBC cost then reduces to the solution of a single n*n 
Riccati equation while determination of the gains requires the solution of an 
additional (n - m) x (n - m) Ricattl equation. This compensator is of general 
interest because it is of sufficient dynamic order to allow arbitrary 
placement of the closed loop poles [7] and is of specific interest here 
because the optimal MODOBC cost represents a lower bound on the optimal SOFC 
cost. (This follows from the fact that any compensator of order zero is a 
degenerate case of a compensator of order n - m, assuming n > m). 

Before proceeding with derivations of other bounds on the optimal SOFC 
cost, it is useful to put this reduced order compensator problem in perspective 
relative to two other reduced order formulations. First, note that no 
explicit control penalty is Included in the cost functional considered here. 

As noted in [1], no such direct penalty is required because as the compensator 
gains increase, so does the contribution of the observation noise to the state 
covariance. In particular, lack of an explicit control penalty does not 
result in a singular problem in this case as it does in the standard linear 
regulator problem. If an explicit control penalty is included in the cost, 
the minimum variance compensator formulation of Martin and Bryson [9,10] is 
obtained. Specifically, the cost functional considered in [9,10] is 

J = lim E{x'(t) Q. x(t) + u’(t) Q„ u(t)} (7-11) 

t-x» 1 c 

providing the designer greater freedom to control the magnitude of the compen¬ 
sator gains at the expense of specifying additional design parameters. The 
principal disadvantage of this compensator is that it becomes ill-posed in 
the constant gain limit since the cost becomes infinite if the control is 
allowed to contain direct feedthrough from a plant output contaminated by 
white noise. Furthermore, since the SOFC problem considered here does not 
penalize the control, the resulting optimal cost represents the best 
performance achievable with any constant gain compensator. Also, since the 
optimal gains are finite, they provide a physically useful indication of the 
difficulty of achieving such performance. 

The other reduced order compensator formulation with which the present 
formulation should be compared is the Levine-Athans-Johnson formulation 
[4,5,6]. The compensator structure considered in this case is the same, but 
the plant driving noise and observation noise are replaced by initial state 
uncertainties and the cost is the same as that for the linear regulator 


87 



problem. Consequently, the resulting optimal compensator depends on the 
Initial state statistics which may be more difficult to specify than 
observation and driving noise statistics. Also, at least partially as a 
result of the greater number of design parameters involved, the necessary 
conditions for this compensator are more difficult to solve than the SOFC 
necessary conditions. 

7.3 A Lower Bound on the Optimal SOFC Cost 

As noted in the last section, the optimal MODOBC cost is a lower 
bound on the optimal SOFC cost. However, it is possible to derive a more 
easily computed lower bound by noting that the necessary conditions can be 
solved in the special case where B is a nonsingular nxn matrix. This 
approach also has the advantage of providing additional Insight into the 
nature of the SOFC problem and, in fact, motivates the upper bounds 
developed in the next section. Specifically, if B - exists and P - ^ exists, 
determination of the optimal SOFC gain reduces to the solution of 


M = -B -1 X C’W -1 (7-12) 

o o 


AX + X A' + V 
o o 


X c’w -1 cx 

o o 


0 


(7-13) 


Note that in this case, the optimal gain M q is independent of the cost 

weighting matrix Q provided only that it is sufficiently well-behaved for the 
solution of equation (7-9) to be positive definite. This would appear to 
imply that the effect of the optimal gain M is to minimize all of the 
asymptotic state covariances independently. 

An even more significant observation is that the asymptotic state 
covariance X q is independent of B itself and only depends on the existence of 

B - *. Because of this, for any given matrix Q the same optimal cost J q is 
achieved by applying the optimal gain (which does depend on B), around 

any open loop plant described by the same A and C matrices and subject to the 

same disturbances V and W, provided only that B exists. In addition, if 
Jg(M) represents the cost associated with applying the gain M around the open 

loop system (A, B, C, V, W), then for any B and § 

BM - BM -> Jg(M) - Jg(M) (7-14) 


A. 

provided A, C, V, W and Q are the same in both cases. Now suppose M* is the 
optimal SOFC gain for the plant (A, B, C, V, W) and that B is not necessarily 
invertible, but that B is. Thus, we can define an "equivalent gain" MT for 



the plant (A, B, C, V, W) by 


M = B BM 


(7-15) 


so that 


BM + = BM => J S (M + ) = Jc(M*) 

D D 


(7-16) 


However, since B * exists, the optimal gain M* associated with (A, B, C, V, W) 
can be computed from (7-12) with the optimal cost computed from (7-13) as 

J q = tr[QX Q ) (7-17) 

Clearly, then, we must have 

J q < Jg(M + ) = Jg(M*) (7-18) 


so that J q represents a lower bound on the optimal SOFC cost for the plant 
(A, B, C, V, W) for any actuator influence matrix B. 

Note that if B 1 does not exist, the computation of the lower bound 
cost J q is easier than the computation of the MODOBC cost since it is not 

necessary to first transform the plant into standard controllable form. If 

B * does exist, it follows that m = n so the MODOBC problem and the SOFC 
problem are identical. In this case, the problem is effectively in standard 

controllable form (i.e., the required transformation is z^(t) = B ^x(0> so 
no additional work is necessary. Finally, note that the covariance equation 
(7-13) is the same algebraic Riccati equation that is used to compute the 
Kalman filter gains required to estimate the state of the open loop plant. 
Consequently, many theoretical and numerical results are available pertaining 
to its solution. For example, it is known [11] that if the open loop plant 
is stabilizable and detectable with respect to the noise sources, there 
exists a unique positive semidefinite solution X q that is positive definite 

if and only if the system is completely controllable from the plant driving 
noise. 


7.4 Suboptimal SOFC Compensators 

The simplifications of the optimal SOFC problem resulting when B * 
exists suggest that it may be profitable to cast the problem in terms of a 

generalized inverse [12,13] of B when B * does not exist. To do this, note 
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that the gain expression may be written as 


M = -B + XC'W _1 


(7-19) 


where the generalized inverse B is given by 


B + = (B'PB) _1 B'P 


(7-20) 


and is in fact a weighted least squares inverse [13] of B. That is, defining 
S as the positive definite square root of P (assuming P is positive definite),, 
i.e., 


S'S = P , S > 0 


(7-21) 


the vector x^ = is the minimum S-norm approximate solution of the 

inconsistent equation Bjc = Consequently, an alternate way of viewing the 
optimal SOFC design problem is as that of minimizing the error 


e - | |Bu + XC’W'^I |g * | |SBu + SXC’W" 1 ^ j (7-22) 

subject to the constraints 

u = (7-23) 

(A + BMC)X + X(A + BMC)’ + V + BMWM'B’ = 0 (7-24) 

P(A + BMC) + (A + BMC)'P + Q = 0 (7-25) 


where S'S « P. Physically, equation (7-23) just specifies the form of the 
control law, equation (7-24) defines X as the closed loop state covariance 
and equation (7-25) relates the norm j|*|L t0 t * ie cost weighting matrix Q. 

Note also that if B exists, it is imnediately apparent that the minimum 
error e = 0 is achievable independent of the norm ||*||g- This formulation 

of the problem thus makes immediately clear the independence of the optimal 
gains from the design parameters noted in Section 7.3 when B is non¬ 
singular . 

This formulation of the optimal SOFC design problem also suggests 
several suboptimal control problems that may be useful in obtaining either 
upper bounds on the optimal SOFC cost or initial gain estimates for 
iteratively computing the optimal SOFC gains. For example, since the 
principal difficulty in computing the optimal gains is the coupling between 






equations (7-24) and (7-25), one possible approach would be to assume a 
least squares weighting S and minimize the error defined in equation (7-22) 
subject to constraints (7-23) and (7-24). Defining P = S'S, the solution 
to this problem is 

M l = -(B'PB) -1 B , PX^C , W _1 (7-26) 

where the closed loop covariance satisfies the equation 

AX X + X X A' + V - X 1 C'W -1 CX 1 = -(I - B(B'PB) -1 B , P)X 1 C'W -1 CX 1 (I - PB(B'PB) _1 B') 

(7-27) 

A particularly convenient choice of weighting matrix S would be the n*n 
identity matrix, reducing the problem to a standard least squares minimi¬ 
zation. 

Another possible approach to decoupling equations (7-24) and (7-25) 
would be to assume a value X for the closed loop state covariance and minimize 
the error defined in equation (7-22) subject to constraints (7-23) and 
(7-25). The solution to this problem is 

M 2 = -(B'P 2 B) _ X B'P 2 XC'W~ 1 (7-28) 

where P 2 satisfies the equation 

p 2 a + a'p 2 + q - p 2 b(b , p 2 b) _1 b'p 2 xc , w _1 c - c , w -1 cxp 2 b(b , p 2 b)' 1 b’p 2 * 0 

(7-29) 

with the least squares weighting matrix S 2 defined as the positive definite 
square root of Pj. A particularly convenient choice of covariance matrix 
would be the lower bound covariance X q defined by equation (7-13) since the 
suboptimal solution then becomes exact when B is nonsingular. 

Finally, a third suboptimal SOFC problem results if both the least 
squares weighting matrix S and the closed loop covariance matrix X are fixed. 
In particular, assuming S - S and X ■ X, the suboptimal gain is given by 

M 3 - -(B , PB) _1 B'PXC , W~ 1 (7-30) 
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A A A /V 

where P * S S as before. Note that if S is taken as the nxn identity matrix 
this approach ia very similar to Kosut's minimum norm suboptimal compensator 
for the Levine-Athans-Johnson reduced order compensator [14]. 

7.5 Conclusions 


Motivated primarily by the exact solution of the optimal SOFC problem 
when the actuator influence matrix B is nonsingular, bounds have been derived 
for the optimal SOFC cost for plants with singular B matrices. These bounds 
satisfy the inequality 

min (J ,J ) S max (J ,J ) < J < JT < J < J ,., (7-31) 

O — O T TT TTT 

where J_ is the optimal MODOBC cost for the plant (A, B, C, V, W), J q is the 

optimal SOFC cost for the same plant but with a nonsingular B matrix, J is 
the optimal SOFC cost and J , J , J ,,, are the costs associated with the 

i T"T “TTT 

suboptimal gains M^, M 2 , and M^ ranked in order of increasing cost. It 

should be noted that if B is nonsingular, it follows immediately that 

J = J = J = J, and that if X is taken as the solution X of equation (7-13) 
o—l o 

in computing M 2 and M^, and are also equal to the optimal cost. 

In terms of computational difficulty, the cost J q is clearly the 

easiest to evaluate, requiring the solution of a single Riccati equation. If 
the plant is in standard controllable form, the optimal MODOBC cost J_ is 
equally easy to compute, but if it is not (as is normally the case), the 
evaluation of J_ can be considerably more difficult. Of the upper bounds, 
the cost J is the easiest to evaluate, requiring the solution of a single 

J A 

Riccati equation if X is available. Evaluation of is somewhat more 

difficult since it requires the solution of equation (7-27) which is not a 
standard Riccati equation. Similarly, the cost should be still more 

difficult to evaluate because of the more complex nonlinearity of 
equation (7-29). Note, however, that all of these bounds on the optimal 
SOFC cost are more easily evaluated, assuming they exist, than the optimal 
cost itself, since they involve the solution of uncoupled equations. 

In terms of applications, three possibilities are apparent. First, 
the lower bound J q and the upper bound should be useful in preliminary 

evaluation of optimal SOFC performance for particular applications. This 
would be particularly convenient if the lower bound covariance X q were used 

in the computation of the suboptimal cost J^. Secondly, an examination of 

the suboptimal compensators M^ and M 2 suggests a possible computational 

algorithm for evaluating the optimal SOFC gains M. Specifically, if 
reasonably efficient algorithms could be developed for solving 
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I 

fl 


equations (7-27) and (7-29), A it might be reasonable to approach the problem 
by first fixing either P or X, solving the appropriate equation to^obtain 
X or P, respectively, and using these to derive updated values of P or X 
from the other equation. Finally, any of the suboptimal gains M^, M 2 or M^ 

could be used as a starting estimate for any iterative calculation of the 
optimal SOFC gains, as has been done with Kosut's suboptimal compensator in 
evaluating the Levine-Athans-Johnson reduced order compensator gains [15]. 
Because of its relative ease of evaluation, M^ appears to be a particularly 
attractive candidate for this application. 


LIST OF REFERENCES 

1. Johnson, T.L., "Optimization of Low Order Compensators for Infinite 
Dimensional Systems", Proceedings of the 9th IFIP Symposium on 
Optimization Techniques , Warsaw, Poland, September 1979. 

2. Johnson, T.L., "Minimum-Variance Fixed-Form Compensation of Linear 
Systems", Proceedings of the 17th Conference on Decision and Control , 
January 1979, San Diego, California. 

3. Interim Technical Report; Actively Controlled Structures Theory , 

Vol. 1, C. S. Draper Laboratory Technical Report R-1249, April 1979. 

4. Levine, W.S. and Athans, M., "On the Determination of the Optimal 
Constant Output Feedback Gains for Linear Multivariable Systems", 

IEEE Trans. Auto. Control , 15, 1, 1970, pp. 44-48. 

5. Johnson, T.L. and Athans, M., "On the Design of Optimal Constrained 
Dynamic Compensators for Linear Constant Systems", IEEE Trans. Auto. 
Control , 15, 6, 1970, pp. 658-660. 

6. Levine, W.S., Johnson, T.L. and Athans, M., "Optimal Limited State 
Variable Feedback Controllers for Linear Systems", IEEE Trans. Auto. 
Control , 16, 6, 1971, pp. 788-793. 

7. Blanvillain, P.J. and Johnson, T.L., "Specific-Optimal Control with a 
Dual Minimal-Order Observer-Based Compensator", Int. J, Control , 28, 

. 1978, pp. 277-294. 

8. Choi, S.S. and Sirisena, J.R., "Computation of Optimal Output Feedback 
Gains for Linear Multivariable Systems", IEEE Trans. Automatic Control , 
Vol. AC-23, p. 760, August 1978. 

9. Martin, G.D. and Bryson, A.E., Jr., "Attitude Control of a Flexible 
Spacecraft", AIAA Guidance and Control Conference, August 1978, 

Palo Alto, California. 




93 






10. Martin, G.D., "On the Control of Flexible Mechanical Systems", 
Stanford University Department Aero/Astro, Report No. 511, May 1978. 

11. Kwakernaak, H. and Sivan, R., Linear Optimal Control Systems , Wiley- 
Interscience, New York, 1972. 

12. Ben-Israel, A. and Greville, T.N.E., Generalized Inverses, Theory and 
Applications , Wiley, New York 1974. 

13. Strang, G., Linear Algebra and Its Applications , Academic Press, 

New York, 1976. 

14. Kosut, R.L., "Suboptimal Control of Linear Time-Invariant Systems 
Subject to Control Structure Constraints", IEEE Trans. Automatic 
Control , Vol. AC-15, pp. 557-563, October 1970. 

15. Petkovski, D.J.B. and Rakic, M., "On the Calculation of Optimum 
Feedback Gains for Output-Constrained Regulators", IEEE Trans. 
Automatic Control , Vol. AC-15, pp. 557-563, October 1970. 


94 



SECTION 8 


CALCULATION OF OPTICAL LINE-OF-SIGHT 


8.1 Introduction 


For the ACOSS test problem //2 we need a line-of-sight model giving 
the law of displacement of the image of a target when the structure deforms. 
The purpose of the model is to give line-of-sight changes of the same order 
of magnitude as those produced by current DARPA large wide-angle optical 
systems. 

We assume that the mirror surfaces are simply displaced and maintain 
their nominal shape. We also neglect any effect of light redistribution in 
the image. 

The structural model given is shown in Figures 8-1 and 8-2. It was 
assumed that the base optical design would be close to an existing DARPA 
design. The mirror shapes are off axis sections of rotationally symmetric 
surfaces which are coaxial. As we want only a first order model and neglect 
the influence of light distribution in the image of a point, the asphericity 
of each mirror was disregarded. 

8.2 Results 


The line-of-sight change is the sum of the changes due to the 
displacement of each of the components of the optical system taken indi¬ 
vidually. For a given mirror or focal plane displacement, the location of 
the intermediate and/or final images are calculated using the well-known 
(see Ref. 1) first order optics equations 


1 




( 8 - 1 ) 


m 



( 8 - 2 ) 


relating the distance d^ between the image and the mirror, the distance 

between the object and the mirror (both distances being taken along the 
optical axis), the radius of curvature r of the mirror and the 
magnification m. The equations (8-1) and (8-2) utilize the vertex of the 
mirror as origin; therefore, when a mirror is displaced, the results are 
given in relation to the displaced positions, and this is taken into 
consideration in the transfer to the next mirror. 
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TERTIARY 



DETECTOR 



Using the paraxial optics equations (8-1) and (8-2), expressions 
relating optical element movements and the resultant longitudinal defocus z^ 
and transverse defocus (x„, y_) are given for the three-mirror ACOSS 

r t 

system. Translations of the mirrors are represented by the displacements 
of the vertices, and rotations are defined for each mirror with the pivot 
point located at its own vertex (on the optical axis). Movements of the 
focal surface are also considered; x^. and y^ are the coordinates of the 

image point, and the pivot point for rotation is again at the intersection 
of the focal plane with the optical axis (see Figure 8-3). Line-of-sight (LOS) 

error is determined by the transverse defocus (x , y ). 

F F 

The pertinent variables are defined as follows: 


Rp = radius of curvature of the primary 

R^_ = radius of curvature of the tertiary 

t^ = axial distance between the primary and the secondary 


t 2 = axial distance between the secondary and the tertiary 

We will use the letter a as a running index which stands for the letter p 
for primary, s for secondary, t for tertiary and F for focal plane, x. 
will represent the displacement of the element a along the ith axis 1 
(see Figure 8-3) and 0 will represent the counter-clockwise rotation 

of the element a about the ith axis (looking towards the origin). With 
these notations the defocus components (x^,, y^,, z^) will be given by 




x. + x. - R 0 + 28- 

lp It p 2p 2s 


+ X lt - X 1F + 9 3F y l 


/R \ I R 

l^ 4 t l/ 2e 2t T 2T - R t 


[ —X— + x. + R 0. - 20, ( -rE- + t, ) + 20, T 1 - 

2p t p lp Is \ 2 1/ It J 2T - R 


(8-3) 


+ X 2t " X 2F ~ 0 3F X I 


(8-4) 


(X 3 P " 2x 3s + X 3t } 


" y l 0 lF + X I 0 2F 


+ X-. - x. 


[R - R + 2(t. + t-)] 
p t 12 


2 3t 3F 
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(8-5) 











where 


T 


R 

_ 2 . 

2 


+ t, 


+ t 


2 


The LOS error is given by 


LOS 



where f is the focal length of the system and is expressed as 


( 8 - 6 ) 


(8-7) 


R R 

t P 

2(2T - R ) 


The magnitude of the LOS error is given by 



(8-8) 


(8-9) 


The location of the image plane is a function of the parameters of the 
three-mirror optical system. The distance BF between the tertiary mirror 
and the image plane (see Figure 8-3) is given by 


BF 


R t T 


2T - R 


( 8 - 10 ) 


and for all cases the focal surface must be placed at this location, or, 
alternately the radius of curvature of the tertiary should be calculated 
using 


2T » BF 
R t 'T + BF 


(8-11) 


so that the final image coincides with the focal plane. 

In the above expressions, all radii of curvature and axial distances 
should have positive values. 

It should also be reiterated that all rotations of the individual 
elements are defined with the pivot points on the optical axis. The justi¬ 
fication for this definition lies in the fact that all three elements of 
the optical system are coaxial, i.e., the centers of curvature of these 
elements lie along a common axis. 
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As a final note, the (x , y , z ) expressions were derived with first 
r r r 

order paraxial optics. As a consequence of this approximation, the angle of 
incidence of the object radiation and all displacements and rotations of 
the individual elements must be small. If this condition is not met, not 
only will higher order terms become appreciable, but the focal length and 
consequently the backfocus BF will change and aberrations of the optical 
system will no longer be negligible. Therefore those equations cannot be 
realistically used in an open attitude control loop. 

8.3 Application to the Test Case #2 

For the test case til corresponding to the structural model given on 
Figures 8-1 and 8-2, the primary and tertiary mirror vertices are coincident 
and the focal plane and tertiary mirror (flat) are coplanar. 

The numerical values of the basic variables which have been used are: 

R 

P 


BF = 

T, R t and f are calculated using Equations (8-6), (8-11), and (8-8). 
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SECTION 9 


DRAPER MODEL #2 


9.1 Introduction 


In order to assess the performance, sensitivity, and hardware 
requirements of the various active structural control methods which are 
being developed, a universal system model is required. This section contains 
a complete description of Draper Model #2 which is a simple but realistic 
evaluation model. The design of this system was driven by the desire to 
incorporate certain attributes into the overall system characteristics. 

The desired features were 

• structural design based on realistic sizes and weights 

• a simple unclassified optical system with associated 
performance measures and tolerances 

• a set of disturbances typical of equipment vibration and 
attitude control (slew). 

The resulting model is described in detail in Section 9.2 along with the 
resulting mode shapes and natural frequencies of the system. Section 9.A 
contains a description of the line of sight (LOS) performance measure 
including theory and implementation. The tolerances on the optical system 
and a set of disturbances are given in Section 9.5. A perturbed model 
which can be used to evaluate the sensitivity of the control system is 
given in Section 9.6. 

9.2 Structural Design 

The example problem is shown in Figure 9-1. It consists of a 
flexible optical support structure and an isolated equipment section which 
contains the solar panels. 

9.2.1 Optical Support Structure 

The optical support structure consists of the upper mirror support 
truss, lower mirror support truss, and the metering truss which maintains 
mirror separation. The finite element model is shown in Figure 9-2. It 
contains 35 node points and 117 beam elements. The structure is configured 
as a truss, but it is assumed that all joints allow a full moment connection. 
Thus, both bending and axial stiffness are included in the model for all 
structural members. All members are assumed to be graphite epoxy hollow 
tubes. The sizes and section properties are given in Table 9-1. The mass 
of this subsystem including the structural mass, mirrors, and equipment has 
been lumped at 18 node points as shown in Figure 9-3. The total mass of 
the optical support structure is 5084 kg. 
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EQUIPMENT SECTION 









model No. 2—Finite element model. 
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UPPER SUPPORT TRUSS - NODE LOCATIONS 
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UPPER SUPPORT TRUSS - MEMBER CONNECTIVITIES 
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9.2.2 Equipment Section 

The equipment section consists of a central rigid body and two 
symmetric solar panels. The support structure for the solar panels is a 
hollow graphite epoxy tube. The structural properties of the solar panels 
are not included in the model. The inertia properties are lumped at 6 node 
points on the support structure. Details of the finite element model and 
the inertia properties of the equipment section are given below in 
Figure 9-4. 



Rigid Body Properties _ Solar Panel Properties _ 

Node 44 Nodes 48, 52, 53, 57 Node 50. 55 


m = 3500 kg m = 90 kg m = 180 kg 

I = 2100 kg-m 2 I = 270 kg-m 2 I = 540 kg-m 2 

xx XX XX 

I * 2100 kg-m 2 

yy 6 

I = 4200 kg-m 2 


116 





9.2.3 Isolation System 

The optical support structure and the equipment section are connected 
by a passive isolator at three points acting in all three translational 
directions. The isolator is designed to reduce the amount of force trans¬ 
mitted from the equipment section to the optics. This is accomplished by 
designing an isolator with the natural frequency much lower than expected 
disturbance frequency. In addition, the isolator frequency should be lower 
than the optical structure bending frequency. A passive isolator consists 
of a spring and a viscous dashpot in parallel as shown in Figure 9-5. 


K 



PRECISION SECTION 


EQUIPMENT SECTION 


_1 / K < M 1 + M 2 » " 

211 V M 1 M 2 

Figure 9-5. 

If the isolator frequency is lower than the lowest natural frequencies of 
the two sections, they can be modelled by their inertia properties only. 

For this system, the isolator was sized to have a frequency of 0.5 Hz. The 
dashpot size has been specified as it is assumed to be part of the overall 
control design. The spring constants are given in Table 9-2. 

9.3 Structural Analysis 

A modal analysis was performed using the NASTRAN finite element 
program. The input deck is listed in Table 9-3. This model contains 
84 dynamic degrees of freedom, thus 84 frequencies and mode shapes were 
extracted. Table 9-4 gives a list of frequencies and descriptions of the 
first fifty mode shapes. 

9.4 System Performance 

The performance of this system is measured by the ability of maintain 
line-of-sight (LOS) rotation (jitter) and defocus within specified tolerances. 
These quantities can be computed using simple linear equations with the 
movement of the mirrors as inputs. The motion of the mirrors is a function 
of the movement of the support points. This section presents the equations 
necessary to compute LOS error from support point displacements. 


F. i 
isol 
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ISOLATOR SPRING CONNECTIVITIES 


TABLE 9-2 


NODE A 

DOF 

NODE B 

DOF 

4 

1 

42 

1 

4 

2 

42 

2 

4 

3 

42 

3 

3 

1 

46 

1 

3 

2 

46 

2 

3 

3 

46 

3 

6 

1 

47 

1 

6 

2 

47 

2 

6 

3 

47 

3 












Table 9-3 
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Table 9-4. Frequencies and mode shapes 



Isolator X-rotation 






























9.4.1 LOS Error Algorithm 

The equations (9-2), (9-3), and (9-4), relating optical surface 
motion and the resultant LOS rotations about the X and Y axes 
below with notations and orientation consistent with the structural 
analysis. The movement of each mirror is defined to be the displacements 
and rotations in the global X, Y, Z directions at a point on the sphere 
(primary, tertiary) or plane (focal plane) which intersects the optical 


(global 

Z) 

axis 

. The LOS rotation and defocus 

are given by 




LOSX = Y/F 

L0SY = X/F 

Defocus = Z 

(9-1) 

where 

F 

_ 

8.051 





X 

= 

A. t-X + X - 
1 p t 

R*0Y + A *0Y - 

p 2 s 

2T*0Y t ] + X - X f 

(9-2) 


Y 

= 

A. f-Y + Y + 
1 p t 

R*9X - A *0X + 

p 2 s 

2T*0X ] + Y - Y f 

(9-3) 


Z 

= 

A.JZ - 2Z + 
3 p s 

V + - z f 


(9-4) 

where 

A 1 

= 

0.2987 





A 2 

= 

93.90 





A 3 

= 

0.0892 





R 

= 

53.9 





T 

= 

66.95 





The terms X , Y , Z , 9X., 0Y., 9Z. refer to the translations and rotations 

1 l ill i 

in the global X, Y, and Z directions of the primary (p) and secondary (s), 
tertiary (t), and focal plane (f). 

9.4.2 LOS Error Algorithm Implementation 

The motion of each mirror is based on the displacements of the 
support points and its position relative to the optical axis. It is assumed 
that the mirrors are kinematically mounted in order to decouple their dynamic 
motion from the structure. A typical kinematic mount is shown in Figure 9-6. 


128 





p 


Y 



Figure 9-6. Kinematic Mount 


Point A is supported in X, Y, and Z directions. Point B in the Y and Z 
directions, and Point C in the Z direction. Using small displacement theory, 
the displacement and rotation of any point P can be found. 


where 



> 

ii 

Ax a + z p A6 y - 

y P A6 z 

(9-5) 


Ay p - 

Ay A - Zp A9 x - 

V 6 Z 

(9-6) 


Az p = 

Az A + y P A6 x " 

X p A6 y 

(9-7) 

A0 

X 

= (Ax c ' 

- Az a + x c (Az a 

- Az B )/x B )/y C 

(9-8) 

A0 

y 

= (Az a 

- Az b )/x b 


(9-9) 

A0 

z 

= (Ay B 

- Ay A )/x B 


(9-10) 


The locations of the support points and point P for each surface in this 
system are shown in Figure 9-7 and are summarized in Table 9-5. 
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Surface 

Node A 

Node B 

Node C 

X 

P 

n 

Z 

P 

Primary 

34 

35 

30-28 


10.0 

0.0 

Secondary 

40 

- 

- 

1 1 

- 

- 

Tertiary 

27 

29 

33-32 


-3.0 

0.0 

F.P. 

11 

9 

40 

-4.0 

-5.0 

0.0 


Table 9-5 


It should be noted that since node C of the primary and tertiary is not 
located at a structural node point, the displacement can be found by inter¬ 
polating between the two nodes given in the table. Also, all six degrees of 
freedom of the secondary correspond to the six degrees of freedom at node 40. 
The location of point P is in the local coordinate system of each mirror 
which are all parallel to the global coordinate system. Using equations 
(9-5) - (9-10) and the information in Table 5-5, a displacement transformation 
matrix (DTM) which relates the motion of each mirror to the displacements of 
the support points in the global coordinate system. 


If modal analysis is employed to compute the displacements, the LOS 
error and defocus for each modal vector can be computed beforehand and the 
total errors found by modal superposition. 

✓ LOSX | 

* LOSY I = [4> lqs ] (u) (9-11) 

where l^Qg] Is a matrix which relates LOS error and defocus to modal 

amplitude. It is computed by first computing the mirror motions due to each 
mode shape using the DTM and then applying equations (9-1) through (9-3). 

T 

For Draper model #2 f^Qg] is given in Table 9-6 for modes 1-50. It was 
calculated using mode shapes which were normalized to unit generalized mass. 

9.5 Parameter Variations 


In order to assess the sensitivity of the control system to changes in 
the natural frequencies and mode shapes, a set of perturbed frequencies and 
mode shapes can be used. This perturbed model can be obtained by either 
varying the frequencies and mode shapes directly or by changing selected 
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Table 9-6 


MOOE 

LOSX 

IOSY 

DEFOCUS 

IOSMAG 

FREQ 

**** 

««*# 

**** 

******* 

****** 

«**» 

1 

0.2705D-04 

0.37780-03 

0.93130-09 

0.37870-03 

0.0001 

2 

-.15800-04 

0.84260-06 

-.18630-07 

0.15820-04 

0.0001 

3 

0.49270-03 

-.31330-04 

-.96040-09 

0.49370-03 

0.0001 

4 

0.11660-03 

0.11550-03 

0.43660-10 

0.16410-03 

0.0001 

5 

-.35020-04 

-.81100-03 

-.2721D-08 

0.8117D-03 

0.0001 

6 

-.78350-03 

0.46760-04 

0.69650-09 

0.78490-03 

0.0001 

7 

-.26410-06 

0.32630-03 

-.4904D-07 

0.32630-03 

0.1455 

8 

0.34390-06 

-.21410-03 

-. 1099D-05 

0.21410-03 

0.2632 

9 

-.15050-05 

-.34520-06 

-.5238D-05 

0.15440-05 

0.3173 

10 

-.22930-03 

0.44070-06 

-.1467D-06 

0.2293D-03 

0.3329 

U 

0.60410-05 

0.10660-05 

0.15590-04 

0.6135D-05 

0.4432 

12 

0.73590-03 

-.67680-05 

0.19060-05 

0.7360D-03 

0.5779 

13 

0.23300-05 

0.87310-03 

-.42760-06 

0.87310-03 

0.5814 

14 

-.10270-05 

-.23620-03 

-.4831D-07 

0.23620-03 

1.2238 

15 

0.37590-05 

-.15190-03 

-.7763D-05 

0.1519D-03 

1.3002 

16 

-.15090-03 

-.14130-05 

0.93640-05 

0.15890-03 

1.3475 

17 

0.27380-04 

-.27770-02 

-.4955D-04 

0.2777D-02 

1.7209 

1C 

0.2590D-05 

-.49760-06 

0.2634D-05 

0.26370-05 

1.8187 

19 

0.11510-04 

-.23640-05 

-.54930-06 

0.11750-04 

1.8187 

20 

-.35970-13 

-.35160-14 

0.3030D-14 

0.3615D-13 

1.8892 

21 

0.22370-04 

0.13760-01 

0.90280-05 

0.1376D-01 

2.3635 

22 

0.13450-03 

0.27240-05 

-.15090-04 

0.1346D-03 

2.9895 

23 

-.43550-05 

-.26020-03 

0.25000-05 

0.26020-03 

3.1759 

24 

0.5573D-06 

-.45360-03 

-.425G0-06 

0.4536D-03 

3.3873 

25 

-.18140-13 

-.43980-15 

0.26010-14 

0.18150-13 

5.1617 

26 

-.33950-04 

-.70130-06 

0.47470-05 

0.33960-04 

5.2603 

27 

0.10380-14 

0.72640-16 

9614D-16 

0.10400-14 

7.8769 

28 

-. 4275D-03 

-.17010-02 

-.19500-01 

0.1753D-02 

8.1168 

29 

0.29160-03 

0.47750-06 

-.363CD-05 

0.2916D-03 

8.3600 

30 

0.18680-01 

-.11030-04 

0.10720-02 

0.13680-01 

8.5706 

31 • 

-.38220-04 

-.2919D-06 

-.3755D-05 

0.38220-04 

8.8135 

32 

-.44940-04 

0.7780D-07 

0.9546D-06 

0.44940-04 

8.0135 

33 

0.37620-05 

0.69950-05 

-.52400-06 

0.9750D-05 

11.3462 

34 

-.12250-02 

0.1070D-02 

0.39820-03 

0.2236D-02 

11.4978 

35 

0.10 Q GD-Q1 

-.56620-01 

-.1593D-01 

0.10990-01 

12.7258 

36 

0.20970-04 

0.74750-03 

-.6865D-04 

0.7478D-03 

13.5832 

37 

-.71430-03 

-.10110-01 

0.11420-02 

0.10140-01 

13.7141 

38 

0.39530-02 

0.59350-02 

1084D-02 

0.71310-02 

14.1604 

39 

-.14140-01 

0.32130-02 

-.49660-02 

0.14510-01 

15.6523 

40 

0.19070-02 

0.10440-01 

-.C294D-03 

0.10620-01 

16.0724 

41 

-.67800-02 

-.22200-02 

0.1413D-01 

0.7134D-02 

16.5248 

42 

0.38530-02 

-.95400-03 

-.13750-01 

0.39690-02 

16.7453 

43 

0.11320-03 

-.64050-02 

-.1739D-02 

0.84060-02 

17.1555 

44 

0.32440-02 

-.1166D-02 

0.71740-02 

0.3448D-02 

17.8203 

45 

-.64390-03 

-.14460-01 

-.1545D-02 

0.1447D-01 

19.0713 

46 

-. 2266D-01 

0.1669D-01 

0.52170-01 

0.2614D-01 

23.7716 

47 

0.90530-02 

0.37110-01 

-.22930-01 

0.33200-01 

24.4140 

48 

0.268S0-03 

0.12050-01 

-.3317D-Q2 

0.1206D-01 

25.908 > 

49 

-.14440-01 

-.10420-01 

0.23620-02 

0.17810-01 

26.3625 

50 

0.94160-14 

-.62570-14 

-.14850-14 

0.11310-13 

26.4292 
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physical parameters which will cause a variation in the modal characteristics. 
Since the effects of modelling and manufacturing errors are of primary 
concern, the latter method was chosen. Two sets of physical parameter 
variations are presented here. Modal analyses including the modal LOS 
error were conducted using NASTRAN. The resulting natural frequencies and 
mode shapes are given and compared with the original system. 

Two sets of parameter variations are summarized in Table 9r7. The 
changes are confined to the optical support structure and involve mass 
distribution changes on the upper optical support truss and stiffness 
changes in the metering truss. These parameters were chosen based on their 
impact on LOS error in several critical modes. The amount of variation is 
driven by the desire to vary the natural frequency by 10%. 

Dynamic analyses were conducted on both perturbed models using NASTRAN. 
The resulting natural frequencies are given in Table 9-8 along with those 
of the original model. The modal LOS errors about the X and Y axes and 
the overall magnitude are given for all three cases in Table 9-9. The 
results show that both sets of parameter variations had more impact on the 
higher frequency modes because the upper support truss and the metering truss 
contain a higher percentage of the total energy in these modes than in the 
lower modes. 

The perturbed models presented here are, of course, two point designs 
from an infinite number of possible parameter variations. In order to fully 
evaluate the stability limits of a particular control system design, it may 
be necessary to evaluate several perturbed models. Another type of variation 
which must also be considered are the changes in amplitude and frequency 
of the system disturbances. 

9.6 System Disturbances and Tolerances 

9.6.1 Description 

The disturbance model consists of two sinusoidal forces which are 
applied to a point on the optical support structure and a point on the 
equipment section. The location, direction magnitude and frequency are 
given in Figure 9-8. These forces simulate on-board vibrating equipment. 

It is assumed that both forces will act simultaneously. 

The tolerances for system line of sight error and defocus are 


LOS = 

1.0 x 

10" 6 

radians 

Defocus = 

.50 x 

cn 

l 

o 

meters 
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Table 9-7 


• Mass changes 


Node 

Original 

Mass 

Case It2 - 2A 

Case //4 - 4A 

27 

375.0 

350.0 

375.0 

28 

375.0 

350.0 

350.0 

29 

375.0 

350.0 

350.0 

30 

375.0 

350.0 

350.0 

32 

500.0 

500.0 

550.0 

33 

500.0 

550.0 

450.0 

34 

250.0 

300.0 

300.0 

35 

250.0 

250.0 

200.0 


• Stiffness changes 


Memb. 

' 

NA 

NB 

Section Property If 

Case 2A 

Case 4A 

76 

8 

14 

1 

1 

80 

9 

15 

1 

1 

81 

11 

17 

l 

1 

85 

13 

19 

1 

2 

99 

15 

32 

1 

2 

101 

17 

33 

1 

2 


• Section Properties 

If 1 - 25cm dia. * 0.10 cm round 
tube 

A = 7.823 x 10" 4 m 2 

I = 6.063 x 10 -6 m 4 

J = 1.213 x 10 -5 m 4 
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It 2 - 25cm dia. x 0.025 cm round 
tube 

A = 1.962 x 10 -4 m 2 
I = 1.529 x 10' 6 m 4 
J = 3.059 x 10 -6 m 4 























Table 9-8 







































Table 9-9 





































































































































































NODE 37 DOF 3 
F - 200 sin (20*t) 







9.6.2 Open Loop Performance 

To demonstrate that structural control is required, the open loop 
performance Is shown In Figure 9-9. In this analysis each mode is assumed 
to have 0.1% critical damping and there is no additional isolator damping. 
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SECTION 10 


PRELIMINARY ASSESSMENT OF DRAPER MODEL 


10.1 Introduction 

The Draper model described in Section 9 was analyzed with respect to 
the two given sinusoidal disturbance functions. The steady state line-of- 
slght, LOS, amplitude error is determined and then compared to the LOS 
performance requirements. The amount of modal damping required to meet the 
LOS performance is determined. The preliminary results show that modal 
damping alone is not sufficient to suppress the sinusoidal disturbance 
impact on the LOS error. 


10.2 Steady State LOS Error 

The amplitude of the steady state LOS error is determined by eval¬ 
uating the modal response at t = °°. The modal equation is 

Hi + 2 Wi + Vi - ♦!, d f 

where 

rj^ is the modal amplitude for the ith mode 

is the modal frequency of the ith mode 

is the damping ratio of the ith mode 

<p is the disturbance influence at location 
’ 'd' to the ith mode 

f is the disturbance forcing function 


For this particular case, the disturbance forcing function is assumed 
sinusoidal. 


f = A sin £lt 

Evaluation of the modal response at t * 



where 


6 - 


-1 


( 


t * 0 

yields the steady-state solution 


. - sin (flt + 6) 

« 2 ) 2 
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The steady state amplitude of the modal response is defined as 


AMPL 



(10-1) 


The steady state LOS error amplitude is obtained by multiplying the above, 
AMPL, by the LOS error influence LO g) for the ith mode. 


£ L0S 


*i,L0S ^ 


(10-2) 


10.3 Sinusoidal Disturbance Impact on LOS Error 

As described in Section 9, two sinusoidal disturbance functions were 
used to degrade the LOS 


f^ = A^ sin Jl^t 


f 2 = A 2 sin f^t 


A l = 200 N 
= 10 He 
at node #37 

A 2 = 400 N 
n 2 = 5 Hz 
at node #46 


The performance requirement for the LOS error is 


£ i,LOS 


( 1 y radian 
1 y radian 
500 y meters 


) 


| LOS error about X-axis 
, LOS error about Y-axis 
[ defocus error in Z-direction 


The asterisk in Table 10-1 shows that 19 of 50 modes have a steady-state LOS 
amplitude greater than the performance requirement. Note that modes 1 to 6 
are the rigid body modes. The defocus error always meets the performance. 
The definitions for the headings in Table 10-1 are as follows: 


LOSXj 


LOS error about X-axis due to disturbance f 
j = 1,2 (see equation (10-1,2)) 


j ’ 


LOSYj 


LOS error about Y-axis due to disturbance f. ; 
j = 1,2 (see equation (10-1,2)) ^ 








Table 10-1. Steady-state line-of-slght amplitude as functions of 
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HKWOHKlHWHrtfjrt«JNr'ir)flNO(rW' , 'a'OMO -«IOlflwa-^^0 ; ONp|OHOM < ONO«»nN-TON 

OT-.^N-0NWlftN0C00 , 'CNUI>0l/lMftC0«0<0>tl0l0 s -tWlftN<0lrtinc0NNN^lft»ft«<C^NNN^N$N«5O 
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DEFOCUSj 


LOS. 

J 


LOS defocus error in Z-direction due to 
disturbance f ; j = 1,2 (see equation (10-1,2)) 


~\j (LOSXj ) 2 + (LOSYj ) 2 ; j = 1,2 


LOS 


V(L0S1)' 


+ (LOS2) 2 


10.4 Modal Damping Requirements 

The damping requirements for each mode are determined by solving 
equations (10-1), (10-2) for and specifying the desired performance 
amplitude, e^ 0 <,. 


C. 


_1 _ 

2w.fi 



i,d^i,L0S 

£ L0S 




2 


Table 10-2 shows the results for the disturbance f^ at node #37 while 
Table 10-3 is for disturbance f^ at node #46. Ignoring the rigid body, 

modes 1 to 6, it can be seen that there are numerous damping ratios » 1. 

In order to obtain some insight, a frequency response was determined (Fig. 10-1) 
between disturbance f and the L0SX. Note the solid line at 62.8 rad/sec 
(10 Hz). Since this line intersects the frequency response at a location 
removed from the peaks, just adding damping will only affect the magnitude 
of the response peaks and not the general shape of the frequency response 
in the neighborhood of the disturbance frequency. 

10.5 Summary 

Preliminary results show that the defocus error is not degraded by 
the two sinusoidal disturbances. Ignoring the rigid body, 17 modes degrade 
the LOS error about the X and Y axes. Furthermore, the addition of damping 
to individual modes is not sufficient to suppress a sinusoidal disturbance 
impact on LOS error. 
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Table 10-2. Damping ratio required to meet LOS performance (f^) 









DAMPING RATIO 

K03E 

FREQ 

AMPl 

L05X 

LOSY 

OEFOCUS 

X 

Y 


(r/s) 


ir) 

ir) 

(m) 



1 

0.000382 

-.67020-04 

-.13150-03 

-.25340-07 

-.62470-13 



o 

0.CC04C9 

.52710-03 

-.83280-06 

.4442D-C9 

-.93190-11 



3 

0.000420 

-.59620-04 

-.29320-07 

.18660-08 

.57260-13 



4 

0.000436 

-.41970-04 

-. 43940-03 

-.43480-08 

-.18320-14 



5 

0.000460 

.17750-03 

-.62173-08 

-.14400-06 

-.43310-12 



6 

0.QC0472 

.13493-03 

-. 1C570-C6 

.63090-03 

.*4240-13 



7 

0.9144C6 

-.76720-04 

.20260-10 

-.25030-07 

.37*30-11 



S 

1.653311 

.40860-04 

.14260-10 

-.e75CD-03 

-.44910-10 



9 

1.593596 

-.262*0-03 

.39470-09 

.90560-10 

.13740-03 



10 

2.091*05 

.52410-04 

-.12020-07 

.23100-10 

-.76680-11 



11 

2.734732 

.39970-03 

. 2415D-06 

.42*10-09 

.62300-03 



12 

3.630322 

-.15740-03 

-.11520-06 

.10653-03 

-.30000-09 



13 

3.653243 

-.17993-03 

-.41910-09 

-.15710-06 

.76920-10 



14 

7.639052 

.71230-05 

-.73120-11 

-.16330-08 

-.3441D-12 



15 

6.169607 

-.1*020-04 

- .60240-10 

.24340-08 

. 12s4D-09 



16 

3.466374 

.67323-04 

-.10770-07 

-.95730-10 

.63493-09 



17 

10.S13C30 

.43660-04 

.11960-06 

-.12130-06 

.21630-C3 



18 

11.427216 

.40210-05 

.10410-10 

-.20010-11 

.10590-10 



19 

11.42/431 

-.33220-05 

-.33260-10 

.73560-11 

.18253-11 



CO 

11.669540 

.11603-13 

-.41730-27 

-.40510-23 

.35143-23 



Cl 

14.650258 

-.72330-03 

-.16180-07 

-.994 *0-05 

-.65300-03 


19.77 

C2 

18.733633 

-. 24303-24 

-.32700-03 

-.66210-10 

.26683-09 



C3 

19.955C63 

.18600-04 

-.81020-10 

- .45410-03 

.4*5iD-10 



C4 

21.232741 

.25500-04 

.14433-10 

-.11750-07 

-.11030-10 



C5 

32.431671 

.15933-14 

-.23910-23 

-.7CC7D-30 

.4144D-29 



Co 

33.051529 

.23603-05 

-.97100-10 

-.CCOeD-21 

.13570-10 



27 

49.492233 

-.66520-16 

-.71100-31 

-.49720-32 

.65380-32 



C8 

50.995:69 

.7652D-03 

-.3271D-06 

-.13010-05 

-.14020-06 


.1750 

C9 

52.5C7542 

-.13510-04 

-.39390-03 

-.64503-11 

.49130-10 



30 

53.e:csi5 

-.65880-03 

-.12310-04 

.73010-08 

-.70620-06 

1.899 


31 

55.3766C2 

.74620-C6 

-.26520-10 

-.21723-12 

-.26020-11 



32 

55.376617 

.26410-05 

-.12770-09 

.22110-12 

.27130-11 



33 

71.2*0451 

-.16970-05 

-.63260-11 

-.15270-10 

.£8540-12 



34 

72.242 c v>6 

.40493-03 

-.49600-06 

.75720-06 

.16120-06 



35 

79.958420 

-.69410-03 

-. 7615D-C5 

.39303-06 

.11C6D-04 

1.838 


3 j 

65.345627 

-.46400-04 

-.97330-09 

-.34*90-07 

.31850-08 



37 

66.163513 

.62470-03 

-.44630-06 

-.63170-05 

.71360-06 


2.003 

33 

63.972637 

-.9951D-04 

-.393hD-C6 

-.59050-06 

.107^0-06 



39 

93.346039 

-.60270-04 

.11350-05 

-.25790-06 

.39860-06 

.2990 


40 

100.935977 

-.65340-04 

-.12460-06 

-.6324D-Co 

.14990-07 



41 

103.328430 

.1307D-03 

-.82600-06 

-.29000-06 

.18460-05 



42 

1C5.214C53 

-.24333-03 

-.93750-06 

.23210-06 

.33450-05 



43 

107.791077 

-.20660-03 

-.23350-07 

.17360-05 

.3592D-06 


.8038 

S4 

112.016337 

.17540-03 

.56900-06 

-.2C49D-06 

.12580-05 



45 

119.828273 

.12330-03. 

-.62600-07 

-.16550-05 

-.19S3D-06 


1.080 

46 

149.3-51542 

-.23160-05 

.52493-07 

-.38650-07 

-.1208D-C6 



47 

153.397614 

.20210-04 

.1331D-06 

.75010-06 

-.46350-06 



43 

162.7SG466 

-.23870-04 

-.7762D-08 

-.34300-06 

.95790-07 



49 

165.640626 

-.22533-03 

.32520-05 

.23480-05 

-.53660-06 

3.492 

2.398 

— 50 

166.055479 

.10660-15 

.1CC40-29 

-.66720-30 

-.15840-30 
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Table 10-3. Damping ratio required to meet LOS performance (f 2 ) 


DAMPING RATIO 


MODE 

FREQ 

AMPL 

LOSX 

LOST 

DEFOCUS 

X T 


<r/«) 


(r) 

(r) 

(m) 


1 

0.000332 

-.5536D-03 

-.15110-07 

-.21100-06 

-.52020-12 


Z 

0.000409 

.42300-02 

-.66830-07 

.3S640-C8 

-.7C79D-10 


3 

0.00C420 

-.87630-03 

--4318D-06 

•2746D-07 

. 64170-12 


4 

0.000436 

-.43030-03 

-.50170-07 

-.49700-07 

-.18760-13 


5 

0.000440 

.14450-02 

-.50730-07 

-.11750-05 

39420-11 

.21050+05 

6 

0.000472 

.17140-02 

-.13430-05 

.6016D-07 

.11960-11 

.29330+05 

7 

0.014426 

.25060-02 

-.6618D-09 

.81750-06 

-. 1229D-09 


a 

1.653811 

-.33470-04 

-.13430-10 

.32390-03 

.42290-10 


9 

1.9935 q 6 

-.30250-03 

.45530-09 

.10460-09 

.15670-06 


10 

2.091605 

.57210-03 

-.13120-06 

.25210-09 

-.63910-10 


11 

2.784732 

-.20560-02 

-.12420-07 

-.21920-03 

-.32050-07 


12 

3.630322 

-.175«D-02 

-.12940-05 

.11900-07 

-.33530-0S 

3.503 

13 

3.653243 

-.21730-03 

-.50630-09 

-.13970-06 

.92930-10 


14 

7.6S9C52 

.95470-02 

-.96090-08 

-.22550-05 

-.46120-09 

3.883 

15 

8.169607 

.11330-03 

.42600-09 

-.17210-07 

-.67970-09 


16 

6.466S74 

-.22020-01 

.35000-05 

.31120-07 

-.2062D-06 

5.771 

17 

10.613030 

-.4S153-06 

-.13200-10 

.13360-03 

.23620-10 


16 

11.427218 

-.14400-02 

-.37310-03 

.71680-09 

-.37940-08 


19 

11.427431 

'.476o0-03 

.54330-08 

-.11270-03 

-.26180-09 


20 

11.869540 

-.80460-11 

.28050-24 

.23310-25 

-.24360-25 


Z 1 

14.850258 

-.56650-03 

-.12630-07 

-.77S2D-05 

-.51140-08 

6.348 

22 

18.733630 

.54340-01 

.73110-05 

.14310-06 

-.S2C2D-C6 

3.891 

23 

19.955063 

.11280-03 

-.49110-09 

-.29340-07 

.231^0-09 


24 

21.2S2791 

-.35«30-01 

-.20030-07 

.16300-04 

.15300-07 

6.497 

25 

32.431671 

-.75850-10 

.13760-23 

.3336D-25 

-.19730-24 


26 

33.051529 

-.84530-01 

.287C0-05 

.59290-07 

-.4012D-06 

.1366 

27 

49.492233 

.10640-12 

.11040-27 

.77260-29 

-.10230-28 


23 

50.955C69 

-.45130-05 

.19290-03 

.76753-03 

.6SCCD-09 


29 

52.527542 

.83250-02 

.25740-05 

.42140-03 

-.32100-07 

1.273 

30 

53.850S15 

-.14520-03 

-.27120-05 

.16050-03 

-.1556D-C6 

1.425 

31 

55.376 6 0 2 

.20360-03 

-.79730-08 

-.60890-10 

-.76340-09 


32 

55.376617 

-.59850-03 

.26910-07 

-.4659D-10 

-.57170-09 


33 

71.290451 

-.31230-07 

-.11750-12 

-.26090-12 

.16360-13 


34 

72.242966 

. 794'*0-C5 

-.97330-03 

.14S6D-07 

.31630-03 


35 

79.953420 

-.33340-04 

-.36810-06 

.16990-07 

.53440-06 


36 

85.345627 

.43780-02 

.51820-07 

.32720-05 

-.30050-06 

3.659 

37 

86.168518 

.31040-03 

-.22170-06 

-.31350-05 

.35460-06 

3.538 

33 

68.972637 

.28010-04 

.11030-06 

.16630-06 

-.30380-07 


39 

93.346039 

.12270-04 

-.17360-06 

.39430-07 

-.60950-07 


40 

100.535977 

-.10640-04 

-.2C26D-07 

-.11110-06 

.24400-03 


41 

103.825430 

-.53260-05 

.36110-07 

.1182D-07 

-.75230-07 


,42 

105.214050 

.43300-05 

.16680-07 

-.4130D-03 

-.59520-07 


*43 

107.791077 

.711SD-05 

.60590-09 

-.59330-07 

-.12360-07 


44 

112.018337 

-.4S43D-05 

-.15720-07 

.56620-03 

-.3477D-07 


45 

119.326273 

-.20770-05 

.13370-03 

.30030-07 

.32100-08 


46 

149.361542 

.29780-06 

-.67480-03 

.4970D-08 

.15540-07 


47 

153.397614 

-.11850-05 

-.10740-07 

-.43990-07 

.27160-07 


46 

162.790466 

-.26350-06 

-.72160-10 

3236D-0S 

.8907D-09 


49 

165.640636 

-.50270-06 

.72580-08 

•524CD-03 

-.1197D-03 


50 

166.059479 

-.12o50-16 

-.11910-30 

.7916D-31 

.1879D-31 
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